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Computer  Programs  for  Three-  Dimensional  Cable 
Problems  in  Tethered -Balloon  Applications 


1.  INTRODUCTION 

In  Reference  1,  six  computer  programs  for  use  in  the  solution  of  various 
tethered-balloon  problems  were  developed  and  documented  therein.  Program  No, 
76.  006  in  Reference  1 provides  a means  of  obtaining  the  tension,  elevation  angle, 
and  space-position  of  a single  tether  cable  of  any  size  and  weight  from  the  balloon 
to  the  ground  for  any  balloon  altitude  and  any  two-dimensional  wind  profile.  The 
wind  profile  from  the  ground  up  to  the  balloon  may  include  winds  of  any  magnitude 
but  all  must  lie  in  the  same  azimuth  plane,  however,  they  can  have  ± signs.  It  was 
therefore  desirable  that  a three-dimensional  case  be  developed  for  use  in  a similar 
type  calculator/computer  to  permit  a completely  realistic  entry  of  atmospheric  wind 
conditions  and  a three-dimensional  evaluation  of  the  cable  geometry  and  other 
physical  parameters. 

Program  Nos.  77.  007  and  77.  007P  presented  herein,  are  three-dimensional 
programs.  They  retain  many  of  the  features  of  Program  No.  76.  006  including  the 
option  of  specifying  a cable  (cylinder)  drag  coefficient  or  calculation  of  a variable 
drag  coefficient  based  on  local  Reynolds  Numbers  at  various  altitudes.  Entry  of 
wind  magnitude  and  azimuth  at  up  to  twelve  different  altitudes  is  permitted.  The 
programs  are  longer  than  No.  76,  006  chiefly  due  to  the  complexities  of  a three- 

(Received  for  publication  15  September  1977) 

1.  Wright,  John  B.  (1976)  Conmuter  Programs  for  Tethered- Balloon  System 
Design  and  Performance  Evaluationr  AFGL-TR-76-0195. 
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dimensional  solution,  the  addition  of  a plotting  routine,  and  the  extra  "housekeeping" 
required  in  dealing  with  azimuth  values.  Program  Nos.  77.  007  and  77.  007P  differ 
by  (a)  the  number  of  optional  ways  to  make  repetitive  computational  runs  and  (b)  the 
plotting  of  parameters  in  77.007P.  A third  program.  No.  77.  007B,  allows  one  to 
retain  a cable  length  computed  with  77.  007,  change  the  winds,  and  then  determine 
a new  balloon  altitude  and  cable  geometry. 


2.  EQUIPMENT 

The  program,  in  three  versions,  documented  herein  was  developed  using  a 
Hewlett-Packard  Model  9810A  calculator/ computer.  All  versions  require  nearly 
full  use  of  the  2036  program  steps  (HP  Option  003),  111  storage  registers  (HP 
Option  001),  the  paper  tape  printer  (HP  Option  004),  and  the  MATH  ROM  (HP  No. 
11210A).  Additional  items  required  are; 

Program  No.  77.007  Printer- Alpha  ROM  (HP  No.  11211A), 

Program  No.  77.007B  Printer- Alpha  ROM  (HP  No.  11211A), 

Program  No.  77.  007P  Printer- Plotter  ROM  (HP  No.  11261A) 

and  an  HP  Model  986 2A  Plotter. 

As  with  those  in  Reference  1,  these  programs  were  written  with  the  idea  that  they 
could  be  adopted  to  other  less  capable  machines  and  they  do  not  necessarily  make 
full  use  of  potentials  of  the  9810A. 

Table  1 defines  the  symbols  used  in  the  program  listings.  Unlike  Reference  1, 
the  listings  (Sections  3.  3. 6,  3. 4.  6,  and  3.  5.  6)  are  not  direct  copies  of  the  output 
tape  listing  where  many  mnemonics  can  have  two  meanings  depending  on  opera- 
tional mode.  Instead  the  listings  here  are  program  forms  showing  meaningful 
mnemonics  both  in  and  out  of  the  alpha  numeric  mode.  In  addition,  the  mathematics 
shown  in  display  register  columns  x,  y,  and  z on  the  forms  will  aid  in  understand- 
ing the  operation  of  the  programs. 
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Table  1.  Program  Codes 


A.  Standard  Mode  — Computation 


In  the  definitions  below,  x,  y,  and  z. represent  the  contents  of  display  registers 
X,  y,  and  z respecavely;  a and  b represent  the  contents  of  memory  registers  a and 
b respectively.  The  mnemonics  and  their  respective  functions  are  shown  below. 


Mnemonic 

TT 

b 

a 

y-^ 

X-* 

1/x 

IND 

xO 

2 

X 

Rt 

l 

t 

f 

X 

+ 

CHG  S 
ENT  E 

CLR 

0 through  9 


CNT 


LABEL 


Function 

71  -*  X 
b -♦  X 

a -♦  X 

y-* memory  register  which  follows 
x-* memory  register  which  follows 
1/x  -♦  X 

Used  for  indirect  addressing 
Puts  the  value  in  the  following  memory 
address  into  x 
x2-.x 

x-»y,  y'+z,  Z-+X 
2'*y*  y'^x,  z-*z 
y -»  x,  X -+  y 
X ■*  y,  y**z,  x-*x 
■*x 
y/x  -*  y 
xy  -♦y 
y - X -♦  y 
y + X y 

-X  *+  X 

Used  to  assign  an  exponent  to  a number 
being  entered  into  x 
.Set  to  0,  X,  y,  z,  a,  and  b 
0 through  9 -*  X 

Used  to  place  a decimal  point  in  a number 
being  entered  into  x 

(Continue)  Used  as  a null  operation  in  a 
program.  Used  in  running  of  a program. 
Used  in  conjunction  with  a following  symbol 
to  Indicate  a position  in  the  program 
memory. 
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PNT 


x<y 


X = y 


X >y 


Go  To 

END 

STOP 

x^ 

In 

X 

e 

arc 

sin  X 
cos  X 
tan  X 
TAB  . 4 
TAB  . 9 


Prints  the  value  in  x.  When  multiple 
PNT's  are  used,  lines  are  skipped 
after  x is  printed. 

If  x<y,  branch  to  address  indicated  by 
number  in  next  4 steps;  if  not,  skip  the 
next  four  steps. 

If  X = y,  branch  to  address  indicated  by 
number  in  next  4 steps;  if  not,  skip  the 
next  four  steps. 

If  x>y,  branch  to  address  indicated  by 
number  in  next  4 steps;  if  not,  skip  the 
next  four  steps. 

Go  to  the  memory  location  specified  in 
the  next  steps, 

l/'sed  as  the  last  step  in  a program; 

Set  point  of  operation  to  Step  0000, 

Causes  program  to  stop  and  permit  entries. 

x^  -*  X 

Natural  logarithm  of  x x 
o x 

Used  in  conjunction  with  sin,  cos  and  tan  keys 
to  obtain  inverse  trigonometric  functions, 
sin  X X 
cos  X ■*  X 
tan  X X 

Common  logarithm  of  x -*  x 

Rounds  number  in  y to  the  power  of  10  indicated 
by  integer  value  of  number  in  x.  Rounded 
number  •*  x,  y unchanged. 
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B,  Alpha  Numeric  Mode  — Printing 

The  alpb  i numeric  mode  is  entered  by  using  FMT  step  twice  and  is  exited  by 

using  FMT  once. 

A through  Z 

A through  Z respectively  is  printed 

0 through  9 

0 through  9 respectively  is  printed 

It 

# is  printed 

-f  or  / 

/ is  printed 

CLR 

This  causes  a carriage  return  or  move  to 

next  line  to  be  printed. 

CNT 

This  causes  a blank  space  in  printing. 

C.  Plotter  Commands 

The  following  combinations  of  key  strokes  serve  to  operate  the  plotter. 

FMT,  t 

Lifts  pen 

FMT,  1 

Drops  pen 

FMT.  1 , t 

Lifts  pen,  scales  coordinates,  moves 

to  coordinates. 

FMT.  1 . ^ 

Scales  coordinates,  moves  to  coordinates. 

drops  pen. 

FMT.  1 . 1 

Symbol  scale  factor  from  x 

FMT.  1 . 2 

Scales  X-coordinate  from  x and  y 

FMT.  1 . 3 

Scales  Y-coordinate  from  x and  y 

FMT.  1 . 4 

Drops  pen,  draws  + at  point  plotted. 

FMT.  1 . 5 

Draws  X-axis. 

FMT.  1 . 6 

Draws  Y-axis 

FMT,  1 , FMT 

Initiates  plotter  alpha  mode 

FMT 

Terminates  plotter  alpha  mode. 

D.  Special  Commands 

FMT,  Go  To 

Automatically  loads  program  card(s) 

FMT,  x-» 

Records  data  in  storage  onto  card(s) 

FMT,  xO 

Automatically  loads  data  card(s) 

To  POLar 

Converts  rectangular  coordinates  in  x and  y 

registers  to  polar  coord,  with  result  placirig 
angle  in  y and  radius  in  x. 
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3.  PROGRAMS  FOR  TETHERED-BALLOON  CABLE,  THREE-DIMENSIONAL 

CASE,  VARIABLE  WIND  PROFILE,  OPTIONAL  INTERNAL  VARIABLE 

DRAG  COEFFICIENT 

3.1  General  Description 

Typical  concerns  in  the  design  of  a tethered -balloon  flight  system  are  the 
ability;  (a)  to  lift  the  weight  of  the  cable;  (b)  to  maintain  the  cable  tension  below 
its  working  limit;  (c)  to  retain  a reasonable  cable  tension  at  the  winch  without  the 
cable  laying  on  the  ground;  and  (d)  to  keep  the  balloon  within  an  acceptable  area 
overhead  under  widely  varying  wind  conditions. 

All  of  the  buoyancy  and  aerodynamic  forces  introduced  into  a cable  system  by 
the  balloon  can  be  summed  and  defined  by  a single  force  and  its  angle.  This  total 
force,  F,p,  and  the  angle,  6 , can  be  computed  by  use  of  either  Program  Nos. 

76.  003,  76.  004,  or  76.  005  in  Reference  1.  These  two  parameters  are  then  treated 
as  inputs  into  either  the  two-dimensional  cable  Program  No.  76.  006,  Reference  1, 
or  these  three-dimensional  programs. 

The  basic  forces  acting  on  the  cable,  in  addition  to  the  total  balloon  force  acting 
at  the  top  of  the  cable,  are  the  cable  weight,  the  aerodynamic  drag,  and  the  ’’e- 
sultant  restraining  force  at  the  ground  winch.  The  weight  of  the  cable  per  ..lousand 
feet  is  specified  while  the  force  at  the  winch  is  part  of  the  problem  solution. 

The  drag  of  the  cable  is  complex  since  it  is  a variable  function  of  the  atmos- 
pheric density,  wind  velocity,  and  cable  diameter.  Since  the  programs  were  in- 
tended for  use  with  balloon  altitudes  of  up  to  66,  000  ft  MSL,  the  effect  of  Reynolds 
Number  on  drag  coefficient  could  not  be  ignored.  Reynolds  Number  is  directly 
proportional  to  cable  diameter,  atmospheric  density,  and  wind  speed,  and  inversely 
proportional  to  atmospheric  coefficient  of  viscosity.  It  was  further  assumed  that 
the  cable  cross-section  would  be  circular  so  that,  in  effect,  the  cable  can  be  con- 
sidered to  be  a cylinder— or  a series  of  cylinders.  Accordingly,  the  program  was 
designed  to  permit  the  user  to  specify  either  a fixed  cylinder  (cable)  held  con- 
stant throughout  the  altitude  range  or  a program  computing  Cp  which  varies  with 
altitude,  wind  velocity,  cable  diameter,  Reynolds  Number,  etc. 

In  concept,  the  cable  is  broken  into  rigid  elements  of  a specified  length,  K,  and 
the  forces  acting  on  this  length  evaluated  to  a net  magnitude  and  angle  with  which 
the  next  lower  element  must  align  and  provide  equal  restraint.  Figure  lA.  Thereby 
a series  of  outputs  is  provided  at  each  of  these  many  points  proceeding  downward 
from  the  balloon  to  and  including  the  surface.  Some  of  these  outputs  are  the  cable 
tension,  space-position,  elevation  and  azimuth  angles,  and  length. 

Because  of  the  three-dimensional  capability  of  the  programs,  not  only  must  the 
wind  velocity  be  calculated  at  each  of  the  elements  but  its  direction  must  similarly 
be  evaluated.  A table  of  altitude-wind  speed-wind  azimuth,  part  of  the  initial  user 
entries,  is  utilized  by  making  straight-line  interpolations  for  each  of  the  two  wind 
parameters  at  each  of  the  element  altitudes.  (See  Section  3.  2.  9. ) 
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During  the  downward  progression  of  calculations,  the  cable  tension  and  pitch 
angle  are  monitored  for  the  condition  of  zero  tension  or  horizontal  cable.  Under 
such  conditions,  the  balloon  has  not  provided  sufficient  lifting  force  for  the  size 
and  weight  of  the  cable  and  the  cable  is  said  to  be  unable  to  reach  the  ground.  For 
a given  balloon,  cable,  and  atmospheric  condition,  a lower  flight  altitude  is  there- 
fore suggested.  If  such  an  event  occurs  just  at  the  surface,  the  cable  is  lying  on 
the  ground.  Hauling  in  some  cable  would  bring  the  balloon  to  a lower  altitude  and 
lift  the  cable  off  the  ground. 

When  the  cable  reaches  down  to  the  surface  where  the  winch  would  be  located  — 
or  where  the  tension  becomes  zero  or  the  cable  horizontal— sufficient  details  are 
presented  as  a final  output  to  allow  construction  of  the  three-dimensional  geometry 
of  the  cable  as  well  as  the  compass  azimuths  of  several  components. 

Table  2 is  useful  as  an  aid  in  selecting  which  of  the  three  programs  is  best 
used  for  a particular  problem.  Programs  77,  007  and  77.  007P,  nearly  Identical 
programs  except  for  ending  rerun  options  and  a plotting  routine  in  the  latter,  are 
based  on  solving  a cable  program  where  the  balloon  altitude  is  fixed  at  a specified 
level.  Both  will  give  identical  answers  to  a given  problem..  Both  provide  the  op- 
tion of  restarting  the  same  program  or  of  lowering  the  balloon  to  other  fixed  alti- 
tudes specified  by  the  table  of  winds  which  was  originally  entered. 

Program  77. 007  offers  two  additional  options;  (a)  a rerun  with  a different 
cable  without  having  to  re-enter  a table  of  winds,  and  (b)  the  ability  to  hold  the 
cable  length  just  calculated  constant,  change  the  w’inds  and  find  where  the  balloon 
may  reach  an  equilibrium  altitude.  This  latter  option  requires  use  of  Program 
No.  77.  007  B.  While  making  use  of  the  same  mathematics  as  the  other  two  programs, 
its  logic  and  details  are  somewhat  different.  The  two  programs,  77.  007  and  77.  007B, 
are  designed  to  be  used  together  and  therefore  each  can  be  used  to  call  in  the  other. 

Program  No,  77, 007  P makes  the  following  plots  on  a single  piece  of  paper: 

(a)  Altitude  vs  H (Y)  Displacement, 

(b)  Altitude  vs  I (X)  Displacement, 

(c)  H vs  I (X  vs  Y), 

(d)  Altitude  vs  Tension, 

(e)  Altitude  vs  Cable  Elev,  Angle, 

(f)  Altltiide  vs  Effective  Dynamic  Pressure.  ' 

As  with  all  programs  so  far  introduced  in  this  series,  these  three  programs 

treat  the  static  condition  and  do  not  attempt  to  consider  the  dynamics  of  balloon 
or  cable  motion. 
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TubU)  2.  HluhliHhtH  of  Pri>|]runiH 


Program  No,  77.007 

Program  No,  77  , 007  P 

Roquirea  2 Curds -4  Sides 

Hrqulres  2 Cards -4  Sldns 

Plus  a Curds -a  Sides 

IJtttn  Cards  with  Density 
and  Cp  Coniitunts  for 

User  ICntriesi  Qalloori  altitude,  Surface  altitude,  Cable  C|j  or  Internal  Compvitation  of  C,-,, 
Cable  diameter.  Cable  weight,  Element  length  to  be  used  In  soUdlon,  aulioon 
total  force  and  angle,  Winds  from  b.aUoonto  surface— magnitvuie  and 
direction  at  up  to  twelve  altitudes 

Program  Cable  'I'enaion  and  Anglo  along  cable  and  at  the  winch.  Cable  Length,  Space 

Solves  for  , PoaUlona  of  points  along  the  cable  and  Anglos  of  the  elements,  Relative 
and  Prints!  positions  of  the  balloon  and  winch,  Sighting  Angles  and  Slant  Range  of  the 

balloon  from  the  winch. 

■ 1 

No  Plotting 

Plots  B Curves —Cable  Tension,  Cable 

Elevation  Angle,  Effective  Dynamic 
Pressure,  X,  and  Y vs  Altitude 
and  X vs  Y, 

Rerun  Options  after  each  Problem  Solution 

0 - New  Prob.  Start  over,  all  entries 

required. 

1 - New  Prob.  but  only  cable  parameters, 

balloon  total  force,  and  angle 
entries— Alt,  / Winds  held  from 
previous  problem. 

2 - Repeat  runs  with  same  cable  with 

balloon  automatically  lowered  to 
each  altitude  in  original  wind  profile. 
Bln  force  and  angle  entered  at  each 
lower  altitude. 

Rerun  Options  after  each  Problem  Solution 

0 - New  Prob,  Start  over,  all  entries 
required, 

2 - Repeat  runs  with  same  cable  with 
balloon  automatically  lowered  to  each 
altitude  in  original  wind  profile. 

Bln  force  and  angle  entered  at  each 
lower  altitude. 

3“  oaFle”len^^^found^n~orl^  T)^.~ 

0 or  1 runs,  change  winds  — ? Now 
balloon  altitude.  Requires  Prog,  No. 
77.007B.  Auto.  Read-In  Call  for 
77.007B  Cards  when  this  Opt. 
specified. 

u 


Program  No.  77.007B 

Requires  2 Cards— 4 Sides 

Data  Cards  not  required 

User  Entries!  Estimated  Balloon  total  force 
and  angle  at  an  estimated  alt. 
to  be  found  by  program  and 
the  new  wind  profile. 

New  balloon  altitude.  Wind 
magnitude  and  direction  at 
that  altitude,  plus  same 
parameters  as  77.007  for 
this  new  condition 

No  Plotting 

Rerun  Options  after  each  Problem  Solution 
0 - New  Prob.  Start  over.  Requires 
Prog.  No.  77.007.  Auto.  Read-In 
Call  for  77.007  Cards  when  this  Opt. 
specified. 

3 - Continue  to  hold  cable  length. 

Change  winds  again. 


Program 
Solves  for 
and  Prints; 
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'i.'i  l>4<volu|)iiiPiil  ol'  l')<|uati()iiii  iiiui  l*r(igi'ainii 

The  basic  fornuilae  ami  description  hereunder  apply  to  all  three  programs. 
However,  there  are  some  details  which  apply  to  77. 007  and  77.  007P  or  are  at  least 
better  understood  by  following  the  logic  In  these  programs.  Any  differences  in 
approach  with  77 . 007B  are  explained  in  Section  3.4.  1. 

These  tether-cable  programs  will  handle  the  three-dimensional  case  where  the 
aisliruith  and  velocity  of  the  wind  may  vary  with  altitude  or  over  the  whole  cable 
length.  Thus  the  balloon,  cable,  and  winch  will  not  necessarily  be  located  in  one 
vertical  plane.  The  wind  velocity  and  direction  at  up  to  12  altitude  points  will  be 
made  a part  of  the  input  data  required.  Straight-line  Interpolations  of  both  velocity 
and  azimuth  angle  will  be  made  part  of  the  program  computations  for  all  inter- 
mediate altitudes.  A wind  vector  will  be  calculated  for  each  incremental  cable 
element  by  determining  the  altitude  of  the  bottom -end  point  of  each  element,  deter- 
mining the  velocity  and  azimuth  of  the  wind  at  that  altitude,  and  assuming  that  these 
conditions  are  constant  over  the  complete  element  length. 

The  objective  of  the  computation  is  to  (starting  at  the  top  of  the  tether-cable 
where  it  is  attached  to  a balloon's  confluence -point  with  a known  balloon  force  vector) 
evaluate  the  cable  tension,  elevation  and  azimuth  angles,  space  position,  etc. . 
moving  downward  to  the  earth's  surface  where  the  cable  terminates  onto  a winch. 

The  cable  can  be  considered  to  be  made  up  of  a series  of  short  and  rigid  cylindrical 
elements  of  length,  K,  attached  by  freely  pivoting  connectors,  Figure  lA.  Each 
element  will  lie  at  an  angle  in  line  with  the  tension  .vector  solved  for  the  element 
immediately  above  it, 

3.2.1  TENSION  VECTOR 

The  first  cable  element  below  the  balloon.  Figure  lA,  is  contained  in  the  ver- 
tical plane  defined  by  the  wind  vector  at  the  balloon  since  it  also  contains  the  balloon 
total  force,  F,j,,  at  an  elevation  angle,  $.  The  wind  vector  Vj,  acting  on  the  first 
element  will  not  He  in  the  same  vertical  plane  and  is  shown  (looking  down  from  the 
balloon)  rotated  clockwise  from  the  balloon  wind  vector  and  initial  vertical  plane  by 
an  angle,  a.  Solution  of  a free-body  diagram  using  only  the  component  of  the 
wind,  the  element  weight,  and  the  tension,  F.^,,  would  provide  a solution  whereby 
the  next  element  below  would  lie  in  the  same  initial  vertical  plane.  That  is  essen- 
tially how  the  two-dimensional  case— Program  No.  7 6.  006— is  handled  with  a single 
vertical  plane  containing  all  elements  of  the  cable  from  the  balloon  to  the  ground 
winch. 

However,  when  the  side  component  of  the  wind,  Vg,  is  included,  the  resulting 
side  force  (drag)  will  rotate  the  bottom-end  tension  vector  out  of  the  initial  vertical 
plane.  A change  in  elevation  angle  will  also  occur  as  in  the  two-dimensional  case. 
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Figure  lA  is  expanded  in  Figure  IB  to  illustrate  all  of  the  forces,  angles,  and 
linear  dimensions.  It  also  indicates  how  the  total  wind  vector,  Vj,  must  be  broken 
into  components  in  order  to  obtain  the  total  aerodynamic  drag  and  then  resolve  it  in 
turn  into  three  manageable  components.  The  wind  vector  can  first  be  divided  into 
two  components;  one  normal  to  the  cable  Vj^,  which  will  be  considered  the  total  drag 
producer,  and  one  parallel  to  the  cable,  V^,  which  will  be  assumed  to  produce 
negligible  skin-friction  drag. 

The  value  of  will  be  used  to  calculate  Reynolds  Number  in  order  to  select 
the  drag  coefficient  from  a stored  table  of  cylinder  Cj^  - R for  the  option  in  which 
the  program  calculates  a drag  coefficient  for  each  element.  The  Cq,  whether  com- 
puted in  this  fashion  or  entered  as  a constant,  is  then  used  with  atmospheric  den- 
2 

sity  and  Vj^  to  calculate  the  total  drag,  D,p,  of  the  element. 

The  total  drag  can  be  divided  into  three  components; 

(a)  in  the  vertical  direction, 

(b)  Djj  one  of  the  two  horizontal  components  lying  in  the  vertical 
plane  containing  the  cable  element, 

(c)  Dg  the  other  horizontal  component  perpendicular  to  the 
aforementioned  vertical  plane. 

The  angle,  0j  is  both  the  angle  of  the  wind  vector,  above  the  horizontal  and  the 
angle  of  the  total  drag,  D,j,,  below  the  horizontal. 

Equations  can  therefore  be  developed  as  follows— the  numerical  subscripts  are 
deleted  at  this  point  for  clarity; 

Given  as  known  quantities:  V,  6,  and  a; 


= V cos  a 

cos  0 = V cos  a cos  9 

„ 2 .,,2  2 ,,2  ,,2  2 2 a 

Vj^  = V — = V — V cos  a cos  9 

(1  - cos^  a cos^  9) 

Vj^  = \^{1  - cos^  a cos^  9)  = V cos  0 

j L 2 2 „ 

cos  0 = il  - cos  a cos  9 

Sin  0=  V^/V  = cos  a cos  0 . 

Considering  the  drag  vectors: 

sin  0 = D^/D,j,— = D,j,  sin  0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
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Therefore 

cos  a cos  9 (9) 

cos  0=  Dj^/D^-H^  = D^  cos  0 (10) 

cos  a = D^/ Dq— ^ Dg  = Dq  cos  a . (11) 

Therefore 

Dpj  = cos  0 cos  a (12) 

sina=  Dg/D^  ^ Dg  = Dj^  sin  « . (13) 

Therefore 

Dg  = cos  0 sin  a . (14) 


Eq.  (5)  is  used  to  obtain  velocity  for  the  Reynolds  Number /C^  extraction  and  in 
the  solution  for  total  drag  in; 

Dt  = -A  (A  = Diam  in  ft  X K)  (15) 

Eqs,  (9),  (12),  and  (14)  are  then  used  to  resolve  D,j,  into  three  components  that  may 
be  introduced  into  the  free-body  diagram  as  follows: 

Summing  the  forces  (aerodynamic,  weight,  tension): 

(a)  Horizontal  Direction  in  Vertical  Cable— Element  Plane, 

F,j,  cos  9+  ~ "^1  ^1  ^1 

(b)  Horizontal  Direction  Perpendicular  to  the  Vertical  Cable— Elemeni  Plane, 

Dg  = Tj  cos  0^  sin  (17) 

(c)  Vertical  Direction, 

F,p  sin  0 = W + D^  + Tj  sin  0^  . (18) 

Eqs.  (16),  (17),  and  (18)  contain  three  unknowns:  Tj,  0j,  and  jS^. 
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From  Eq,  ( 17 ) 

T = 

1 cos  sin  /3j  * 

Substituting  in  Eq.  (16) 

Ft  6 + Dy  = cos  Sj  'si'n  j3j,  ®1  ^1 

Ft 

^1  ■ P “cos  0 4-  u 

1 H 


From  Eq.  (18) 

F^  cos  6 + D„ 

^ 1 ' cos  cos 

Substituting  in  Eq.  (18) 

F_  cos  6 + D„ 

F.,  .in  8 - W - D„  . J.  ^ » 

cos  6 + Dfj 
= cos 

F„  sin  0 - W - D 
01  = arc  tan  - cos-pVT^  ^1 

F.J,  cos  0 + Djj 
'^l  cos  cos  /Jjj  ' 


sin  6^ 
tan  0j 


(19) 


(20) 

(21) 


(22) 


(23) 

(24) 


Eqs.  (21),  (23),  and  (24)  are  utilized  for  the  solution  defining  the  tension  vector 
at  the  bottom  end  of  the  cable  element  with  which  the  nejct  element  aligns.  In  sub- 
sequent loops  through  these  equations  a value  of  at  the  top  of  an  element  would 
be  used  instead  of  F,p  used  for  the  first  element.  The  particular  forms  of  equations 
developed  were  chosen  over  other  possible  forms  to  avoid  indeterminate  solutions. 
For  example,  when  /3  = 0-deg,  the  sin  (3=0  and  the  tension  would  be  Infinity  if 
Eq,  (19)  was  incorporated  into  the  program.  It  is  assumed  that  (3  will  never  equal 
90“  particularly  if  small  K values  are  selected. 
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ERAL  Breakdown  of 
LE  INTO  Elements 
Length^  K . 


Details  of  the  Forces,  Angles,  and 
Displacements  Used  in  the  Solution 
of  a Three-Dimensional  Cable  Problem 


A set  of  axes  is  defined  as  shown  in  Figure  IB  which  will  be  the  reference  axes 
during  the  main  body  of  computations.  Since  the  origin  is  at  the  top  of  the  cable  at 
the  confluence  point  of  the  balloon,  they  are  annotated  as  follows:  (1)  Yg  axis 
positive  ahead  of  balloon  in  the  direction  from  which  the  balloon  wind  is  blowing, 
and  (2)  Xg  axis  positive  90®  clockwise  from  the  Yg  axis— looking  down  from  above 
the  balloon.  The  reason  for  this  particular  set  of  axes  will  be  apparent  in  subse- 
quent discussion. 

As  evident  in  Figure  IB,  the  initial  element  lies  in  the  vertical  plane. 

The  solution  of  6y  and  Tj,  not  only  establish  the  conditions  at  the  top  end-point 
of  the  next  element  below  but  also  define  a new  vertical  plane  containing  this  second 
element.  This  second  plane  is  rotated  by  the  angle  j3^  clockwise  from  the  Yg-Z 
plane  for  the  example  shown;  the  figure  basically  shows  the  positive  sign  conven- 
tions for  angles  and  distances.  This  second  vertical  plane  containing  the  second 
element  becomes  the  plane  to  which  its  wind  vector,  V is  referenced  to  deter- 
mine the  associated  relative  wind  angle,  . A repeat  of  all  the  above  computations 
will  then  provide  solutions  for  ^2’  "^2  therefore  conditions  defining  the 

third  element  down. 

While  the  angle  /3  will  usually  be  small,  the  angle  a can  have  any  magnitude. 

In  these  programs,  a value  up  to  360®  is  permitted  since  the  triginometric  output 
of  moat  calculators  will  handle  a full  360°.  In  the  example  shown  in  Figure  IB,  a 
lies  between  0®  and  90°.  In  this  case  the  vertical  component  of  the  total  drag  is 
downward,  Dg  in  a direction  to  cause  a decrease  in  the  cable  elevation  angle,  and 
Dg  in  a direction  to  cause  a clockwise  rotation  of  the  cable.  These  are  all  positive 
in  the  sense  that  they  were  used  as  illustrated  when  writing  the  free-body  equations. 
As  the  angle  a is  increased  to  90°  and  beyond,  the  three  drag  components  change 
directions  and  at  times  some  disappear  as  illustrated  in  Figure  2.  However,  as 
can  be  checked  by  substitution  of  appropriate  trig  functions  in  Eqs.  (9),  (12),  and  (14), 
proper  values  and  signs  result  for  any  value  of  a as  used  in  these  equations. 

3.2.2  CABLE  GEOMETRY 

Returning  to  Figure  IB,  it  is  apparent  that  while  determining  the  progression 
of  cable  tension  and  angles  working  down  the  cable,  that  the  linear  movement  of 
the  cable  must  be  calculated  and  both  angular  and  linear  data  must  be  properly 
summed  as  computations  proceed.  As  stated  previously,  the  first  cable  element 
below  the  balloon  lies  in  the  vertical  Yg-Z  plane.  The  position  of  the  bottom  end- 
point of  that  element  is  defined  by: 

= K sin  e (25) 

where  is  the  vertical  drop 
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hj  = K cos  e (26) 

where  hj  is  the  horizontal  displacement  from  the  top  of  the  element  in  the  Y direction. 
However,  the  next  elements  move  out  of  the  Y„*Z  plane  and  as  the  cable  re- 

O 

spends  to  side  loads,  the  bottom  end  point  of  each  element  is  displaced  in  both  the 
X and  Y directions.  Similarly,  the  vertical  plane  containing  the  cable  element 
moves  to  a different  azimuth  than  the  next  element  above.  Knowledge  of  this  azi- 
muth is  required  with  the  next  wind  azimuth  in  order  to  define  the  relative  wind 
angle,  a , acting  on  the  element. 


jg  = K sin  Oj  . 

(27) 

g2  = K cos  0,  . 

(28) 

hg  = gg  cos  = K cos  Oj  cos  . 

(29) 

^2  ' ^2  = K cos  sin  0^  . 

(30) 

At  the  next  or  third  element,  further  complications  arise  since  the  above  solu- 
tions would  provide  displacements  along  and  perpendicular  to  the  direction. 
Figure  3,  which  is  an  X - Y projection  of  Figure  IB,  will  clarify  the  following: 
While 

hg'  = K cos  $2  cos  ^2  (3  1) 

hg  a K cos  $2  cos  (^2  + (32) 


or 


hj  = K cos  cos  (^j_j  + /3j_2  + . . . + ^j) 


(33) 


and 


H = hj  + h^_j  + hj_2  + . . . hj  = Fh  . (34) 

While 

ig'  a K cos  ^2  /Sg  (35) 

ig  « K cos  $2  sin  (02  + 0j)  (30) 
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■ jrfl«3iV«rr^."'i'7^«'<'?';{(w,"'J^’"‘" 


or 


i. 

i 


K cos  sin  (/3._j  + "^  “ ' • '*'  ■'^1^ 


(37) 


and 


I = ii+  ii_i  + ii.2  + •••  ^2  ” 

jg  = K sin  02  (39) 


or 


jj  = K sin  0j_j  (40) 

J = 3i  + ji-1  + Ji.2  + • • . + Ji  = • (41) 

Eqs.  C’S),  (37),  and  (40)  are  thus  used  to  determine  the  displacement  of  the 
bottom  end-point  of  one  element  from  that  of  the  element  immediately  above.  Eqs. 
(34),  (38),  and  (41)  provide  the  summation  of  these  distances  to  the  balloon  which 
become  part  of  both  the  element  output  and  the  final  output  when  the  surface  is 
reached. 

The  sum  of  the  /3  angles,  as  used  in  Eqs.  (33)  and  (37),  when  added  to  the  azi- 
muth of  the  wind  at  the  balloon  provides  the  azimuth  of  the  vertical  plane  containing 
the  cable  element.  The  relative  wind  angle,  a,  to  which  the  element  is  being  sub- 
jected must  be  determined.  A positive  convention  is  established  whereby  a is 
always  positive  with  values  between  0®  and  360®.  By  definition: 

a ” Wind  AZ  — Element  AZ  , 


As  in  Program  No,  7 6,006,  wind  (and  density)  conditions  are  found  for  the  altitude 
of  the  bottom  end-point  of  the  element  but  are  assumed  to  exist  as  a constant  over 
the  whole  element  length,  K.  Therefore; 

Qf^  o Wind  AZ  at:  (Zg  - jj)  ^ Wind  AZ  at;  Zg  (42) 

Og  = Wind  AZ  at:  [Zg  - (jj  + jg)]  - Wind  AZ  at:  Zg  - j3j  (43) 

Og  = Wind  AZ  at:  [Zg  ' + J2  -*3^^  " 

or 


= Wind  AZ  at:  (Zg  « Jj)  - Wind  AZ  at:  Zg  - (/3j,j  + /Sj.g  + 
Expressed  In  shorthand  symbols  used  in  some  figures,  this  becomes 
Oj  . AZj  - (j3  g + 51/3 ) , 


‘(45) 


(46) 
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3.2.3  REYNOLDS  NUMBER  ~ DRAG  COEFFICIENT 
Reynolds  number  is  computed  by: 


(47) 


If  a constant  drag  coefficient  is  not  specified,  Cj^  is  calculated  by  one  of  two 
methods  that  are  dependent  on  the  value  of  R.  When  R < 1,  the  Stokes  condition  is 
assumed  and 


C^=  (10.9/R)  / (0.87  - LogR)  . (48) 

When  R>  1,  a series  of  straight  lines  are  used  to  approximate  the  variation  of 
cylinder  Cj^  with  R together  with 

^D  ■ Base  Point  ^ ^Base  Point^  ‘ 

The  following  constants  are  included  in  program  storage  for  Cj^  solution  when  R > 1. 
The  Recall  Code  Numbers,  n^j^,  are  explained  in  Section  3.  2.  8. 


Recall  Code 
Number 

R Region 

^Base 

P 

D Base 

^R 

'’cd"'  ^ 

R > 

1 < 9 

1 

12.  5 

- 10,  0 

'^CD  “ ^ 

R < 

900 

9 

2.  98 

- 1.0 

"CD  *3 

R.  < 

4500 

900 

0.  98 

0.  0 

"CD  * 

R < 

9000 

4500 

0.  98 

0.7308 

"CD  “ 3 

R < 

40, 000 

9000 

1.  2 

0.  0 

"CD  =3 

R < 

50, 000 

40,  000 

1.  2 

- 4. 54 

"CD  = 

R < 

250,000 

50,000 

0.76 

0.3434 

"CD  ' 3 

R > 

250,000 

250,  000 

1.  0 

0.0 

3.2.4  COEFFICIENT  OF  VISCOSITY 

The  Coefficient  of  Viscosity,  is  required  for 

the  calculation  of  Reynolds 

number  in  Eq.  (47).  The  values  of  ^ in  the  1962  Standard  Atmosphere  can  be 
defined  by  two  straight  lines  up  to  86,  000  ft  within  the  accuracy  needed  here. 
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From  Z = 0 to  36.  500  ft!  M = 1.  205“®  - 6.  84164"^ ^ Z Ib/ft-sec  (50) 

From  Z = 36,500  to  66.  000  ft:  M = 0.  95528“®  Ib/ft-sec  . (51) 

3.2.5  DENSITY 

The  atmospheric  density,  p,  is  found  as  a function  of  altitude,  Z,  by  use  of  the 
following  equation: 

In  p/pQ  = a^  Z + a^  Z^  (52) 

where  p^  la  the  density  at  sea  level 

a^  = -2. 813606'®  . 

aj  = -1.777  17“^®  . 

These  constants  were  obtained  from  a fit  of  the  1962  Standard  Atmosphere  up 
through  70,000  ft  with  a precision  considered  satisfactory  for  this  particular  appli- 
cation. Additional  refinements  or  use  of  other  atmospheres  more  typical  of  sea- 
sons or  locations  of  a particular  balloon  flight  could  be  easily  adapted  by  a change  of 
the  two  constants. 

3.2.6  TERMINATION 

The  cable  element  lengths,  K,  whether  entered  or  automatically  made  equal  to 
(Zg-Zg)  / 100,  are  summed  in  Storage  Register  No.  024  at  each  loop  in  the  calcu- 
lation of  end  position,  tension,  etc.  When  an  altitude  for  the  bottom  end-point  of  one 
cable  element  is  detected  below  the  surface  altitude,  provision  is  made  to  go  back 
to  the  altitude  of  the  top  end-point  of  that  element,  divide  K by  10— as  well  as  W 
and  A—  and  then  proceed  again  downward  until  a bottom  end-point  goes  below  Zg 
at  which  time  the  final  printouts  occur  for  the  point  just  above  the  surface.  In  effect 
this  process  provides  a vernier  and  a solution  closer  to  an  exact  value  at  Zg  than 
possible  if  K were  left  unadjusted. 

Therefore  when  resetting  for  the  optional  runs  at  lower  balloon  altitudes,  K,  W, 
and  A are  multiplied  by  10  to  reestablish  their  original  values.  A vernier  of  10 
followed  by  another  10  could  be  incorporated  if  an  extremely  precise  surface  alti- 
tude match  were  desired.  Its  need  would  probably  only  exist  for  a zero-wind  condi- 
tion where  the  cable  elevation  angle  is  very  large  and  the  error  greatest. 

Checks  are  made  of  both  tension  and  cable  elevation  angle  for  positive  values 
before  looping  back  in  the  program  and  adding  another  cable  element.  Should  either 
not  be  positive,  an  appropriate  message  is  delivered  and  the  same  final  printouts 
are  provided  as  when  a surface  condition  is  reached. 
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3.  2.  7 KNU  OF  CAULF-WINCMI 


Wlion  tlu)  altitude  of  u bottom  end-poiUt  of  un  element  reuchex  the  xurfiunt,  the 
run  Is  complete  luui  a finnl  printout  Ih  provided  whleli  ropreaenta  condltionx  at  the 
cubic  winch.  .Similarly,  if  the  tension  uhould  become  zero  and/or  the  cable  become 
horizontal,  the  final  printout  would  roprosont  conditions  at  the  cable  end  at  an  alti- 
tude above  the  actual  surface  altitude,  liut  consider  in  this  discussion  that  the  cable 
has  reached  the  surface  winch. 

Figure  4 illustrates  a complete  Idealized  balloon-to-wlnch  cable  plan  view.  In 
working  the  problem  from  the  balloon  downward  in  altitude,  the  sign  of  the  Vj^^-axls 
is  + ahead  of  the  balloon  and  the  sign  of  the  X^-axis  Is  + 90*  clockwise  from  the 
+ Yy-axia.  Because  of  this,  a simple  transformation  of  axes  to  the  which  locution 
allows  direct  application  of  u conventional  sign  designation  to  the  winch  axes, 
and  Y^,  For  example,  the  sum  of  the  h-diatancos,  H,  is  positive  ahead  of  the 
balloon  and  therefore  when  referenced  to  the  winch  position,  this  distance  would  be 
also  positive  In  the  example  shown.  In  all  but  extreme  cases,  H would  always  bo 
positive.  However,  the  distance  1 can  be  either  + or  - depending  on  relative  wind 
but  would  also  bo  directly  referenced  from  the  winch  and  its  conventional  X-Y  sot 
of  axes. 

Therefore,  the  final  output  data  includes  the  H and  I distances  which  can  bn 
plotted  directly  as  Y and  X distances  on  the  conventional  sot  of  winch  axes  to  show 
relative  wlncli-balloon  positions.  The  straight -lino  distance,  L la  also  printed.  By 
use  of  the  equations  shown  in  Figure  4,  the  azimuth  of  the  balloon  from  the  winch, 
AZp,  is  presented.  By  use  of  the  height,  J and  the  horizontal  distance  out,  L,  the 
elevation  angle  and  slant-range  to  the  balloon  are  provided. 

In  addition  to  the  cable  tension,  length,  and  weight,  the  elevation  angle  or  pitch 
of  the  cable  above  the  horizontal  Is  also  printed.  Use  of  the  equation 

AZ^  = r/3  + /3  - 180“  (53) 

provides  the  azimuth  of  the  cable  leaving  the  winch. 

At  this  point,  machine- or  hand-plotted  data  on  tracing  paper  could  bo  overlaid 
on  a map  of  the  area  with  the  v/inch  point  at  its  known  location  and  the  paper  ro- 
tated to  place  the  north  vector  on  true  north  if  the  exact  geographic  balloon  loca- 
tion were  desired.  To  avoid  this  complexity,  a third  set  of  axes  was  established 
with  the  origin  at  the  winch  and  the  winch  axes  rotated  to  place  the  + Y axis  on  true 
north.  As  shown  on  Figure  4,  distances  X and  Y which  can  be  plotted  directly  on 
a map  or  chaid;  are  calculated  ard  are  the  last  two  items  printed  in  the  final  output. 
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3.  2.  a HOUSEKISISPINCI  DETAILS 


a,  A/iimuth  Mivnugomont  tlu'ovigh  0*  to  300* 

It  was  ilooUtod  to  keep  all  spociflod  unulds  within  tho  0*  to  3B0*  tiosignation 
whother  necessary  or  not  in  order  that  no  ambiguity  exist  particularly  on  a printout 
of  data.  A few  examples  of  detailed  handling  foUowsi 

When  calculating  tho  wind  avslmuth  ut  a specific  altitude,  the  wind  data  entered 
at  tho  nearest  altitudes  above  and  below  are  used  to  obtain  the  rate  of  change  of 
aximuth  with  altltvide  (strulght-line).  As  with  the  wind  velocity,  the  pj-ocess  is 
downward  so  that 


d A / 


A/ 


lower 


- A4, 


upper 


~2 


- Ai4., 


lower 


where  n a NW  « tho  point  number  in  wind-fielii  table  starting  with  NW  » l for 
highest  altitude  point.  .Should  one  ualmuth  lie  in  the  <lth  quadrant  and  the  other  be 
in  the  lat  quadrant,  for  example,  an  improper  value  of  delta  assimuth  would  be 
determined  (u/imuth  angles  are  always  assumed  to  lie  within  tho  < 1110*  included 
angle  between  tho  two  lumwn  axlmuths). 


0“ 


n 

n'^n+l 

Z 11,000 

Z 8000 

AZ  3ia* 

60* 

270 

00 

leo- 


d A/  dO"  - 31.1"  -2a.1» 

or  “ TiTWcrnnnro  “ wit” 

A wind  lu'.imuth  for  Ab  0000  ft 
would  improperly  be  foxuvi  in 
the  3nd  quadrant  at  14.1" 
inste.ad  of  25". 


Similarly,  reversing  the  two  known  points  would  lead  to  an  improper  solution 
at  9000  ft  with  40"  instead  of  350*. 

To  handle  this  problem,  the  following  checks  are  Included  in  the  pi'ugram  logici 


Trial  A Aaimuth,  A A/^  re  ^^n+1  ~ ^^'n 

iao»v 

No  Yes 

I 

IsAAZ,  < -180"? 

' I 

No  Yes 

f ' 1 

AAZoAAZ^  A Aa  a AAZ^  + 360"  A A^  o A A2l^  - 380* 

When  tho  final  aximuth  la  found,  an  angle  outside  tho  0"  to  380“  llmlla  might 
bo  indicated.  To  prevent  any  ambiguity,  it  is  convened  by  tho  following  process! 
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la  A/,  > aBO»V 


N|0 

la  Aii  > 0 V 


Yfla 


tl 


A;^  B A^  + 300* 


Yns 

aI 


B A(^' 


A/  B Ai?.  - 300*  . 


A thii'il  pi'oblam  of  tlUa  natui'o  ai'latia  whan  datai'iDining  tha  nn({lai  Ofi  which  ia 
the  ungle  batwa^^n  the  wind  uulmnth  anti  the  H»in\uth  «)f  the  vai'ticat  plana  containing 
the  cable  element.  The  unglci  a,  ta  ulwaya  Uefinecl  in  the  olockwiae  aenae  aa 
poaitivo. 


« « Wind  A'^  • Kloment  A'/, 


01  H A/  - (r/1  + (i  ij) 

where 


^ B Wind  A^  at  uc  NW  b I . 
If  the  following  wore  preaenti 


rt  r.  aoo*  - 320*  b -lao* . 


A poaltlve  convention  and  printout  ia  maintained  by,- 


Trial  Cl,  iK^  A A - {Tp  -I  p j^) 

la  01^  *-  0“? 

No  Yea 

CK  i 0 i 01^  + 300“. 

Thus  in  the  above  example,  ct  ° i;40“. 

b,  Special  Problem  — Cable  \ ertical 

In  the  usual  totliered-balloon  cable  system,  cable  weight  is  sufficiently  largo, 
relative  to  cable  drag,  that  the  cable  "sags"  or  its  elevation  angle  decreases  with 
decreasing  altitude.  A special  case  can  be  found  by  introducing  a large-diameter 
and  very  ilght-welgbt  cable  into  a strong  wind  field  which  causes  the  cable  to 
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P •■li* mil  ..II 'Sf»n»T9'^a^?'Va’'*^t ' 


1 5 ^ HL  un^.-js^in 


tt\o»'«uw(t  ltd  (tU'VHtlon  imftlw  tliruu^l*  HO*  tU  uomti  uUUutltta,  When  tl»ia  linppens,  an 
amhlj<v\U,y  In  the  emiatiun  ntlUaetl  euuaea  u Ninltlan  anil  false  reversal  to  a negative 
tension  an-.l,  if  unoorreoteit,  a reuiUnit  that  the  tension  is  aero,  To  take  earn  of 
tills  speeial  ease, a snUroutlne  (l-AHl'll.,  nl  is  introihioeil  wlileh  may  bo  examinetl  in 
the  flow  (,'hartiK  It  essentially  oorrecls  the  eompateif  anfile  /i  by  aikliim  1110“  to  fi 
when  tint  eable  noes  through  tlui  vertieal.  The  physleul  mounlnu  of  this  oorreatlon 
Is  apparent  In  the  horizontal  plane  (11-1)  where  0 (joes  thronah  00*,  the  projoetion 
of  the  eable  has  a reflex  of  UlO*  ehanae  of  clireution,  (See  Kiauro  30), 
e,  Storaae  and  Iteeall  Codes 
(1)  Draa  Coefficient 

Tlie  elaht  points  doflnlna  the  draa  eoefflciont  of  a cylinder,  when  HeynolcJa 
Number  U areater  tlian  1.  0,  each  consists  of  three  constants  us  described  in  Sec- 
tion ;i,  3,3,  They  are  stored  in  Storaao  Itoaister  No,  077  throuali  100. 

Uecall  of  the  parameters  la  made  by  use  of  Code  Number  called  n^jij,  A 
value  of  1 Is  assianed  fur  the  group  of  parameters  defined  by  the  smallest  It  or 
Itegisters  No.  077,  0711,  and  070  follosvml  by  a valueof  two  for  the  next  group  and  on 
up  to  elaht,  Kach  of  the  parameters  are  e.vtraeted  wlitm  the  proper  H urea  is  found 
by  the  use  of  Indirect  addressina  utiiUina  these  formulae; 

3 n^,|^  + 74  for  Itjj  , 

3 n^,|j  +73  for  C|j|j » 

3 H(j|j  + 70  for  Kj^  , 

(3)  Wind 

The  wind  profile  may  bo  detined  by  as  few  ua  two  altitude  points  or  as  many 
US  13,  Nach  point  consists  of  throe  usor-ontorod  values;  altitude,  wind  velocity, 
and  azimuth,  Kor  both  entry  and  recall,  use  is  made  of  a t.'ode  Number  callocl  NW, 
l>'or  storage,  use  is  made  of  Indirect  udtiressing  and  the  following; 

3 NW  + 30  for  Altitude  , 

3 NW  + 37  for  Wind  Velocity  , 

3 NW  + 30  for  Azimuth  . 

The  value  of  NW  begins  at  1 for  the  first  group  of  entries  or  highest  altitude  and 
increases  by  1 at  each  lower  point  down  to  and  including  the  surface.  Following 
entry  of  all  points  in  the  wind  field  a value  of  NW  >=  1 is  reassigned  since  the  com- 
putation process  starts  at  the  highest  altitude  and  works  downward. 
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Hocull  of  tho  paramotors  In  Program  No,  76.  006  was  made  by  the  same  formula 
listed  above.  In  these  programs  however,  the  optional  rerun  cycles  with  the  balloon 
at  successively  lower  altitudes  listed  In  the  wind  table  required  a modification.  A 
Itepeater  Code  Number,  r,  is  used  herein  for  several  purposes.  A value  of  0 is 
used  during  the  primary  runthrough  with  tho  balloon  at  maximum  altitude  (llun 
No,  1),  If  tho  option  rorvui  cycles  are  called  for,  r is  indexed  + 3 for  each  succes- 
sive run  at  lower  altitudes.  The  printed  Hun  Number  la  found  by: 

Hun  Number  » 1 + J • 

Hocall  of  tho  wind  parameters  is  made  by  indirect  addressing  of  the  following 
forn'iula: 

3 NW  + 36  + r for  Altitude, 

3 NW  + 37  4 1'  for  Wind  Velocity, 

3 NW  4 38  4 r for  Azimuth  , 

3.  2.  9 WIND  PHOKIL13 

The  usual  available  wind  data  consists  of  the  wind  magnitude  in  knots  and  its 
azimutli  (direction  from  which  tho  wind  is  blowing)  in  degrees  from  true  North  at 
various  altitudes  above  the  surface.  For  a low  altitude  balloon,  the  twelve  available 
altitude- wind  storage  groups  in  the  programs  are  more  than  enough  to  handle  the 
typical  amount  of  wind  data  available. 

For  higher  altitudes  where  more  than  twelve  levels  of  wind  information  is 
available, some  editing,  smoothing,  or  averaging  may  be  necessary.  The  magnitude 
and  azimuth  may  each  bo  plotted  and  points  defining  significant  changes  in  each 
parameter  be  used  to  define  significant  altitude  levels.  As  an  example,  in  Figure  5, 
tho  wind  magnitude  and  azimuth  are  plotted  as  points.  Some  points  can  be  usually 
Ignored  as  dubious  or  offering  little  effect  on  the  total  wind  picture.  Significant 
points  defining  the  wind  niagnitude  variation  with  altitude  as  a series  of  six  straight 
linos  are  found  at  the  following  altitudes: 

28.000  13,500 

23.000  6500 

17.000  3000 

0 . 

Significant  points  defining  the  wind  azimuth  variation  with  altitude  as  a series  of 
seven  straight  lines  are  found  at  the  following  altitudes: 

28.000  9000 

25.000  3000 

22,500  0 . 


Figure  5.  Typical  Wind  Information 


Since  four  of  these  are  the  same  altitudes  defining  the  wind  magnitude,  a total 
of  eleven  altitudes  are  sufficient  to  define  the  magnitude  and  azimuth  variation  as- 
shown  by  the  Table  on  Figure  5. 

The  use  of  straight  line  interpolations  between  entered  wind  magnitudes  and 
between  entered  azimuths  for  each  intermediate  altitude  is  based  on  the  simplicity 
of  this  type  of  operation  and,  in  effect,  the  restructuring  of  the  original  averaged 
wind  profile  described  above.  It  does  not  pretend  to  reflect  the  best  meteorological 
technique  for  precise  wind  analysis.  However,  if  reasonable  care  and  judgement  is 
used  in  selecting  the  input  points  defining  the  two  parameters,  the  resulting  output 
cable  position  and  conditions  will  be  within  the  acceptable  engineering  standards 
used  throughout  the  program. 

As  described  in  Sections  4.  1,  A and  B,  two  adjacent  input  wind  points  should 
never  include  azimuth  angles  that  are  exactly  180°  apart.  In  order  to  avoid  an 
ambiguity  as  to  which  way  a wind  field  is  rotating  between  altitudes  in  such  rare 
cases,  altitudes  should  be  selected  so  that  adjacent  azimuth  angles  differ  by  less 
than  180°. 
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3,  2.  10  MAGNETIC  PROGRAM  AND  DATA  CARDS 
A.  Program  Cards 

Program  No.  77.  007  and  No.  77.  007P  each  require  two  HP  No.  9162-0012 
magnetic  program  cards  utilizing  four  sides  each  program.  As  explained  in  their 
operating  instructions,  several  versions  of  each  program  might  be  desirable.  For 
example  the  following  three  sets  of  cards  for  each  progx-am  have  been  found  useful 
in  avoiding  excessively  long  rolls  of  output  data  when  it  is  not  required: 

(1)  Print  no  intermediate  altitude  plus  Print  final  surface  data  for 


data  from  balloon  to  surface 

(2)  Print  only  Z,  H,  1,  6,  and  T 
values  from  balloon  to  surface 

(3)  Print  all  intermediate  data 
from  balloon  to  surface 


cable  condition  at  the  winch. 

Print  final  surface  data  for 
cable  condition  at  the  winch. 

Print  final  surface  data  for 
cable  condition  at  the  winch. 


Program  No.  77.  007B  also  requires  four-sides  of  two  program  cards.  How- 
ever, no  printing  of  intermediate  altitude  data  is  provided  (see  Section  3.4.  1). 

B.  Data  Cards 

In  order  to  conserve  program  steps,  in  the  development  of  the  programs,  twenty- 
six  drag  coefficient  and  density  constants  were  stored  on  magnetic  data  cards.  There- 
fore, two  HP  No.  9162-0012  cards  (three  sides)  are  required  to  place  the  constants 
in  the  registers  chosen.  As  explained  in  the  operating  instructions,  these  cards  are 
loaded  into  the  machine  in  the  process  of  loading  the  program  cards  and  the  con- 
stants will  be  retained  through  any  number  of  problem  solutions  until  the  machine  is 
turned  off.  One  set  of  data  cards  can  be  used  interchangeably  with  Program  No. 

77.  007  or  Program  No.  77,  007P:  Program  No.  77.  007B  receives  the  constants 
through  Program  No.  77.  007. 

The  two-page  listing  of  the  data  program,  which  follows  hereunder,  inserts  the 
constants  into  the  proper  storage  registers  and  then  provides  a printout  of  all  stor- 
age register  numbers  and  contents  for  review.  Register  numbers  000  through  074 
and  101  through  108  should  contain  zeros.  After  reviewing  the  contents  of  register 
numbers  075  through  100  for  accuracy,  the  data  cards  can  be  recorded  by  the  fol- 
lowing key  strokes; 

END 


Alternately,  if  data  cards  are  not  used, the  constants  may  be  loaded  into  storage 
directly  by  key  punch  each  time  the  program  is  loaded.  The  following  step  numbers 
should  be  changed  as  shown  to  prevent  automatic  Insert  Card  light  and  loading  motor 
turn-on. 


Change  contents  to  CNT: 


No.  77,  007 
Step  Nos. 

()040,  1,  and  2 


No.  77.  007P 
Step  Nos, 
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3.3  Program  No.  77.007 

This,  the  basic  version  of  the  three-dimensional  tether-cable  program  provides 
only  printed  outputs.  Therefore  the  MATH  ROM  and  a PRINTER- ALPHA  ROM 
only  need  be  installed  in  the  HP  9810A  Calculator.  If  plots  are  desired.  Program 
No.  77.  007P  should  be  used, 

3.  3.  1 ACCESS  TO  PROGRAM  NO.  77.  007B 

As  explained  in  Section  3.  3.  2 and  under  Section  3.  4,  this  program,  77. 007, 
must  be  used  as  a first  step  in  entering  Program  No.,  77.  007B.  Once  tied-in,  these 
two  programs  can  be  called  back  and  forth  for  solution  of  many  detailed  problems. 

3.  3.  2 RERUN  OPTIONS 

When  a problem  is  solved  and  the  cable  end  conditions  printed,  the  program  is 
designed  to  offer  the  user  a choice  of  methods  with  which  to  proceed  to  the  next 
problem,  A message  is  printed: 

OPT.  ENT  (OPTION,  ENTER) 

followed  by  a choice  of  four  numbers  and  a STOP. 

When  a completely  new  problem  requiring  reentry  of  different  altitudes,  winds, 
and  cable  parameters  is  next  to  be  run,  the  number  0 is  entered.  The  program  will 
clear  all  but  the  storage  registers  loaded  from  the  data  cards  and  cycle  back  to  the 
initial  printing  of  the  program  number  and  title. 

There  are  situations  where  a series  of  problems  involve  only  changes  in  the 
cable  specifications,  the  element  length,  or  the  balloon  total  force  and  its  angle. 

For  this  situation,  where  the  altitude's  and  the  wind  Held  table  remain  unchanged 
from  one  problem  to  another,  the  entry  of  the  number  1 in  the  above  STOP  will  save 
these  parameters  unchanged  from  the  previous  problem.  The  program  will  clear 
only  summation  registers  and  cycle  back  to  the  initial  printing  of  the  program  num- 
ber and  title.  The  STOPS  normally  used  for  entry  of  altitudes  and  the  wind  field 
will  be  by-passed  although  these  saved  parameters  will  be  printed  at  their  proper 
locations. 

The  normal  runs  (printed  RUN#1  MAX  ALT)  provide  a cable  solution;  (a)  for  the 
balloon  at  a specific  altitude  subjected  to  the  winds  specified  in  the  first  group  of 
wind  entries,  and  (b)  for  the  cable  subjected  to  a variety  of  wind  conditions  specified 
at  the  lower  altitudes  down  to  the  surface.  One  practical  and  sometimes  limiting 
problem  is  concerned  with  raising  or  lowering  the  balloon  through  the  same 
specified  wind  field.  Rather  than  a reentry  of  all  variables  required  at  each  of  many 
lower  altitudes,  a method  is  provided  to  simplify  a series  of  runs  at  decreasing 
altitudes. 

If  the  number  2 is  entered  at  the  above  STOP,  the  program  will  clear  the  sum- 
mation registers  and  retain  the  surface  altitude,  cable  specifications,  and  wind 
field,  and  setup  RUN//2  with  the  balloon  at  the  second  lower  altitude  in  the  wind 
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IttiSU.  ftetiBUHtk  ihB  httUooM  Uitttl  rtml  nuyln  Ht  thin  n«w  tiUUvule  will  dll'ifei' 
fHMti  Un'  ii'MvUnirt  bullmiM  alllUHbt  '*Hl,V  *•«'?  .I/XU'  la  nttetlatj  for  tUBPe  two  entrios, 

A full  iiolulvtMt  U (hoo  iM'vwHloit  (or  UilM  uootUfloo.  Af  itM  oitdi  Iho  projiruin  loopn  lo 
HUNfxX  Aiul  tltH  l)lkUoon  »t  Uio  tiurd  lowor  tillltgdn  lu  Iho  wind  fiold  'lublo,  Th^so 

ttvlomrtlio  loopii*  ronllnvitt  onlll  snrfuoe  Isi  renehtnl,  At  that  palm  tli«  Of*T.  WNU 
id  pi'lntvd  hot  a rlialoo  uCoiili'  0 or  I Id  than  p^r^nltteoi 
Whan  aondliU'viit^  tta  use  of  tnia  Utwarlni;  altUuda  ojielo  option,  It  Is  neuesswi'y 
to  InitUUy  anlai'  a wind  field  twl,'la  wnloh  sp^-3til‘.'H  a sii/fioiant  mimbov  m'  ttltltudos 
to  .mNivldo  0 covnpletfi  oiittli’sls,  AKil\iiieH  hovlna  poidlonlarlji  strontj 

wiinia,  of  aonrovn  in  hnlloon  parformwioo,  would  of  tani'so  ha  invludml  in  tho  Inltiftl 
doflamon  oX  tho  wind  field  a,(tin/i  on  the  oabUn  .Vdinllurl'/,  It  lo  >'ot,'t»asnr,v  to  have 
fho  values  of 'ilie  balloon  total  foinjo,  K^,,,  and  Its  angle,  U,  BvaUob,e  for  entry  wl 
eaob  aUttude  aw  IntUoaled  In  the  INIHIT  OA'J'A  bGllM,  Tivo  use  of  Program  No, 

Yd,  OQd  or  No,  Yft,  00ft,  flefewme  },  for  eweb  of  the  altitudes  Ir,  the  wind  (laid  table 
will  provide  tbo  v,‘,.\ues  o^  (''rj,  and  d, 

'I'l  a effr(,'t  of  a rbange  In  wind  on  the  balloon  oable  systrin  with  IV  juat- 
oabnuated  oabls  lengt,,  bald  flvwd  is  « praetioal  pr'oblem,  'this  qucsilou  van  be 
aiiawai  ed  by  use  of  the  number  X entered  at  the  abovu  OPT  h'NT  S'lXlP  — unly  after 
UA  urlgi.  ,d,  OPT  (!,  or  OPT  I solutloii,  As  tfiv  messafto  indicates,  this  requires 
the  use  of  f'rnKram  No,  77,  007  U whose  eards  should  be  rnady  for  loading  when 
seleovlng  Ihla  option,  U Ibis  r5'pluo.9in0Ul  by  77,  00'/ H is  maik'i  roadmlssloii  to 
77,007  oan  be  waally  made  according  to  the  directions  oontalord  under  77,00711  (see 
heotlon  vS,*},  Pats  card*  are  not  ,\Q''jded  in  these  Interchanges  as  the  cr»?»stant8  in- 
volved are  in'oteotod  by  pixjrain  logic, 
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— , jxVtJ 


u t»  •»? 


1 


\ NQi104« 
u H*40000  >n — :r-v 


lOi,  » ft) 


PNVl 

‘S*- 

•V 

Q<)K 

» Dili  B 

PNVl 

ft, 

a ooa  1 

PNVl 

“s" 

e,, 

nin 

a eoa  6 

h" 

A?« 

' ’•*'  6i 

Oo  \'o  uMteu 

" » 

l.OOR  COVI'i 

No. i0«8 


PNVl 


I’HTl  ej»  AVO  tAtt 


PyAln  9 ’ W 

0 ♦ D„ 


P_  «10»  9 ♦ 0„ 

PNVl  Vv»  “ 

DO 


la  R>JS0000? 


IND  3ll^j,  ♦ 7A 

"o 

INO  3Dj,p  + 7S 

Sb 

INU  3llj.^j  + 76 

<•0  ■'=00  ^ 

-LOG  Rjj) 

(sub  D>  Baud) 

/""Xvoa 

I. 


);b+(i„  -360 


Ift+Bn  - CBtB„ + 


D„»  C„qA 


JmSrfi 


S,  !1, 4 OJMQHATlNfl  INSTUUCTIONN  ANU  NO'IMSS 


KJiiV  STIiOKWa 

fiiuD 

LitiO 


ntioH 


KNTHUiiti 


■■*°(Numi)«i'  of  tloolmiil  i^tnuas  timaU'iad) 
(to  I of  Pi*ogi‘Am  Cttixia 


Continue  inaortlitd  tsUdoH  S,  3,  aud  4 of  Pi'ogi'am  Cat'da 
iKttOi't  Bt(to  I of  Untu  Cai'da 


(Following  houdiug  la  now  printed) 

PROG  H 77.  007 
3 -DIM.  TETHER 
TE.ST 

(X)  (Y)  (Z) 

At  STOP  0-0-0,  Enteri  Test  No/I.  D.  

Following  printing  of  the  test  number  or  1.  D. , Side  1 
of  the  Data  Card  will  load. 

Continue  inserting  Sides  2 and  3 of  Data  Cards 

Note  that  loading  of  Program  and  Data  Cards  will  only  bo 
required  one  time  until  machine  is  turned  off. 

(Following  is  now  printed) 

RUN  #1  MAX  ALT 


At  STOP  1-1-1,  Enter: 


CONT 


JT  or  Cq  * Zg  (ft,  MSL)  Zg  (ft.  MSL) 

♦ Enter  7i  to  use  built-in  cylinder  Cp  variation 
or  enter  Cp  value  which  will  be  held  constant 
throughout  program  run. 


At  STOP  2-2-2,  Enter: 

0 or  K* 

(Element 

Length) 

Wt/ 1000  ft  (lbs)  Diam  (in.) 
(Cable)  (Cable) 

♦Enter  0 to  set  K = — — or  enter  K in  ft 

1 CONT 

At  STOP  3-3-3,  Enter; 

6 (deg) 

(Angle  of  Fy 
to  Horizon)^ 

Fjj.  (lb)  

(Balloon  Tot. 

Force) 

1 CONT 1 

WINDS  (Printed) 

At  STOP  4-4-4,  Enter;  Azimuth  of 

wind  (deg) 

This  entry  must  be  for  conditions 

CONT  1 

Wind  (knots)  Zg  (ft  MSL) 

at  balloon  altitude. 

At  STOP  4-4-4,  Enter: 

Azimuth 

Wind  Z 

Stop  4's  will  repeat  until  the  last  set  of  entries  for  the  surface  condition 
at  Zg  are  inserted.  A total  of  12  sets  of  wind  entries  may  be  made 
including  the  balloon  and  surface  conditions 

CONT  I 
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Tlift  alH)v«  fltttH  of  entries  will  be  prlntofl  oui  In  gi'oupN  whon  CIONT  Ih  uti'uok  ufter 
Sto|)N  It  St  3,  und  4,  Tho  oi'  oomputmi  valuo  will  be  pilnted  in  the  oatie 

ofK. 


Thfl  pi'Ojjram  now  begina  ooinputution  aturting  »t  the  bc^Uoon  and  woi'ka  downward 
one  element  at  a time  to  the  aurface,  The  following  twenty-throe  parametors  are 
printed  for  the  condition  at  tho  bottom  ond-polnt  of  oaoh  elementt  NQTEi  To  avoid 
the  printout  of  any  or  all  of  the  twont,y“throo  parameters  shown  below,  replace  PNT 
with  CNT  at  tho  associated  program  stop  numbers. 


Stop  No. 

0488 

Altitude,  ft  MSL 

0480 

j 

Vert.  Distance,  top  to  bottom  of  element,  ft 

0613 

h 

Horiz.  Oistanoe.  parallel  to  balloon  axis,  ft 

0524 

i 

Horiz.  Distance,  perpendicular  to  balloon  axis,  ft 

0632 

J 

Total  Vert.  Distance,  balloon  to  bottom  of  element,  ft 

0538 

H 

Total  Horiz.  Dist. , Eh,  balloon  to  bottom  of  element,  ft 

0540/1 

I 

Total  Horiz,  Dist. , Ei,  balloon  to  bottom  of  element,  ft 

0545 

ZK 

Total  Element  (Cable)  Length,  ft 

0555 

EW 

Total  Element  (Cable)  Weight,  lb 

0656 

Wd 

Wind  Velocity  at  bottom  of  element,  knots 

0755 

AZ 

Wind  Azimuth  at  bottom  of  element,  deg 

0757 

^ B 

Azimuth  of  Vert.  Plane  Containing  the  Element,  deg 

0778 

a 

Wind  Incidence  Angle  on  the  Element,  deg 

0902 

R 

Reynolds  Number 

1084 

Cd 

Drag  Coefficient 

1098 

q 

Dynamic  Pressure,  Ib/ft^ 

1108 

Dt 

Total  Element  Drag,  lb 

1147 

°w 

Vertical  Drag  Component,  lb 

1159 

Horiz.  Drag  Comp,  in  Vertical  Plane  of  the  element,  lb 

1173/4 

Horiz.  Drag  Comp,  to  Vertical  Plane  of  the  element,  lb 

1201 

^1 

Horiz.  Rotation  of  Tension  Vector  at  bottom  end  or 
Horiz.  Rotation  of  the  next  element's  vertical  plane, 
deg. 

1235 

1293 

^1 

Pitch  Angle  Downward  of  Tension  Vector  at  bottom  end 
or  Elevation  angle  of  next  element  above  the  horizon, 
lb 

1248/9 

1303/4 

^1 

Tension  at  bottom  end  or  at  top  end  of  next  element,  lb 
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Gi'ovips  of  tho  vttluoN  ot'  the  ubove  23  piirunuitora  will  continue  to  bo  printed  for 
points  down  the  civble  vmtll  one  of  the  following  conditions  is  ancounlorodi 

(a)  The  cable  reuchos  the  earth's  surface  ut  Zg—tho  winch  location, 

(b)  The  tension  becomes  itoro, 

(c)  Tho  cable  becomes  horl'/ontul. 

(In  («),  tho  computational  techniques  used  do  not  yield  u precise  Zg  condition.  The 
final  Z will  bo  higher  than  Zg  by  an  amount  less  than  K sin  0/10,  usually  no  more 
than  a few  feet). 

Tho  final  printout  includes  the  abbreviated  names  and  values  of  the  following  param- 
eters. They  describe  Uie  conditions  ut  tho  winch  if  condition  (a)  is  attained  or  at 
the  cable  lower -end  wlrich  is  above  the  surface  if  conditions  (b)  or  (c)  are  indicated. 


(a)  ON  SURFACE  or  (b)  TENSION.  0 or  (c)  CABLE  HOR 


ALT 

HT 

TENSION 

C.ELEV.ANG 

LENGTH 

WT 

V.DRAG 
H.  l.  L 


AZ.  TO  BLN 
B EL.  ANG 
S.R 

CABLE  AZ 
X +E 
Y +N 


Z Altitude,  ft 

J Vertical  Height,  ft 

T Cable  Tension,  lb 

0 Elevation  Angle  of  Cable  above  Horizon,  deg 

TK  Cable  Length,  ft 

m Cable  Weight,  lb 

Total  Vertical  Drag  Component,  lb 
H Tot.  Horiz.  Distance  along  Yg  or  Y^  axis,  ft 

1 Tot.  Horiz.  Distance  along  Xg  or  X^  axis,  ft 

L Min.  Direct  Horizontal  Dist.  to  Balloon,  ft 

AZg  Azimuth  Angle  to  Balloon,  deg 

€ Elevation  Angle  to  Balloon,  deg 

SR  Slant  Range  to  Balloon,  ft 
AZ^  Azimuth  Angle  of  Cable  (Out  of  Winch),  deg 
X X Coordinate  to  Balloon,  ft 

Y Y Coordinate  to  Balloon,  ft 


At  this  point  the  initial  problem  entered  is  solved  with  printing  completed.  If  this 
was  an  initial  run  (RUN  #1  MAX.  ALT),  the  following  is  printed  and  STOP  provided 


to  permit  3 optional  ways  to  rerun  the  program. 

OPT.  ENT 

0 - NEW  PROB 

1 - NEW-SAME  ALTS  / 

WINDS 

2 - LOW  CYCLES 

3 - HOLD  C.  LENGTH 

CHG.  WINDS*  HAVE 
77.  007 B CARDS 
READY 


At  STOP,  Enter:  0,  1,  2,  or  3 in  (X) 
ICONT 
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If  0 is  Entered;  New  Problem— Use  for  completely  new  problem.  All  but  the 
permanent  storage  registers  containing  density  and  drag  coefficient  constants  will 
be  cleared,  program  returns  to  start  with  reprint  of  Number  and  Title. 

If  1 is  Entered:  New  Problem— but  Same  Altitudes  and  Winds— Use  when  only  the 
cable  diameter,  cable  weight,  element  length,  or  balloon  force  and  angle  is  to  be 
changed  In  next  problem.  Summary  storage  registers  will  be  cleared  and  program 
returns  to  start  with  reprint  of  Number  and  Title.  Program  proceeds  on  with  only 
STOPS  2 amd  3 for  associated  entries.  STOPS  1 and  4 are  not  activated  but  the 
altitude  and  wind  field  data  will  be  automatically  printed  at  the  proper  locations 
from  storage  registers  loaded  from  the  previous  problem. 

If  2 is  Entered;  Lower  Altitude  Cycles  — Use  when  analysis  is  desired  with  the 
balloon  lowered  to  each  of  the  altitudes  specified  in  the  wind  field  table.  Summary 
storage  registers  are  cleared  and  the  program  returns  to  start  with  reprint  of 
Number  and  Title.  Instead  of  RUN  #1  MAX.  ALT  the  following  is  printed: 

RUN  # 

2 

B.ALT/AZ 

^B 

Program  then  goes  to  STOP  3 for  entry  of  balloon  total  force  and  angle  at  this  new 
altitude  Zg  which  is  the  second  altitude  originally  entered  in  the  wind  table.  The 
program  will  then  make  a complete  solution  to  the  surface  (or  to  zero  tension  or 
horizontal  cable)  for  this  new  balloon  altitude.  However  at  the  end  of  RUN  #2,  no 
option  is  provided  since  the  program  cycles  on  to  the  third  altitude  point  in  the  wind 
field  table,  sets  up  RUN  #3  and  goes  to  STOP  3 again  for  entry  of  the  two  balloon 
parameters  at  this  nev/  altitude.  Thus  a solution  is  provided  for  each  altitude  in 
the  wind  field  table  except  for  Zg.  When  Zg  is  detected,  the  program  will  terminate 
to  the  OPT.  ENT  printout.  In  this  case  however,  a choice  of  only  0 or  1 will  be 
offered.  If  1 is  chosen  the  program  will  recycle  back  to  the  initial  maximum  alti- 
tude and  properly  retain  the  complete  original  wind  table. 

If  3 is  Entered;  Hold  Cable  Length  and  Change  Winds— Use  this  to  find  balloon  alti- 
tude with  the  cable  length  determined  in  a MAX.  ALT  run  held  constant  and  the 
system  reaction  to  a different  wind  profile  to  be  entered  in  Program  No.  77.  007 B. 
After  the  3 is  entered.  Side  1 of  Program  No.  77.  007B  card  should  be  inserted. 
When  CONT  is  pressed,  the  reading  of  the  4 card  sides  will  then  proceed.  In  this 
case  the  instructions  in  Section  3.4.4  should  then  be  followed. 
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Note  No.  1 - If  incorrect  data  is  believed  to  have  been  entei'ed  do  not  press 
STOP  END  to  restart  prograna.  For  correct  and  safe  clearing 
of  registers  press  following: 

(stop" 

Igo  toI 

m 

E 

m 

After  OPT  ENT  message  is  printed,  enter  0,  then: 

ICONT  ' 

Note  No.  2 - Proof  of  solutions  can  be  made  by  summation  of  the  vertical  forces. 
Summary  of  the  horizontal  forces  is  not  possible  because 
summations  along  two  fixed  horizontal  axes  could  not  be  handled 
within  the  available  apace. 

Therefore:  F,p  sin  6 = T^  sin  + SW  + TD^ 


Note  No.  3 - An  ambiguity  situation  is  set  up  when  a no-wind  condition  at  all 
altitudes  is  introduced,  that  is,  0 = 90®  at  Zg  and  wind  entries 
are  all  0 knots.  In  such  a case; 

(a)  The  slant  range  becomes  infinity  instead  of  being  equal  to 
the  height  or  cable  length  and, 

(b)  The  terms  AZg  and  AZ^  compute  as  180®  opposite  the 
azimuths  used  in  the  wind-field  entries.  Since  the  cable 
elevation  angle  at  tlie  winch  is  90®  for  this  condition,  these 
latter  terms  have  no  meaning— the  balloon  is  directly  over  the 
winch. 
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3.  3.  5 INPUT  DATA  FORM 

Test  No. : Date: Notes: 


INPUT 

77,  007  1 

STOP  NO, 

ITEM 

VALUE 

MAX.  BALLOON  ALTITUDE 

Ft.  MSL 

1-1-1 

SURFACE  ALTITUDE 

Ft,  MSL 

CABLE  Cg 

or  ir  for  Internal  Comp. 

CABLE  DIAMETER 

D 

Inches 

2-2-2 

CABLE  WEIGHT  per  1000  ft 

Lb. 

CABLE  ELEMENT  LENGTH,  KORO 
forK=  (Zg  - Zg)/100 

Ft. 

3-3-3 

BALLOON  TOTAL  FORCE 

ANGLE  OF  TOTAL  FORCE 

F 

T 

e 

■H 

4-4-4 

WIND  PROFILE 

For  Stop  3 
Opt.  Lower 
Alt.  Cycles 

For  Stop  3 

Opt,  Lower 

Alt.  Cycles 

!.■  Zg 

Ft.  MSL 

7.  Z 

Ft.  MSL 

^T 

Lb. 

Wind 

Knots 

Wind 

Knots 

e 

Deg. 

AZ 

Deg, 

AZ 

Deg. 

2.  Z 

Ft.  MSL 

^T 

8.  Z 

Ft.  MSL 

Ft 

^9 

Wind 

Knots 

0 

Wind 

Knots 

6 

AZ 

Deg. 

AZ 

Deg. 

m 

3.  Z 

Ft.  MSL 

Ft 

.9.  Z 

Ft,  MSL 

Ft 

Wind 

Knots 

$ 

Wind 

Knots 

e 

AZ 

Deg. 

AZ 

Deg, 

g 

4.  Z 

Ft.  MSL 

^T 

Lb. 

10.  Z 

Ft,  MSL 

J’t 

Lb. 

Wind 

Knots 

e 

Deg. 

Wind 

Knots 

e 

Deg. 

AZ 

Deg 

AZ 

Deg. 

5.  Z 

Ft.  MSL 

Ft 

11.  Z 

^T 

Lb. 

Wind 

Knots 

6 

Wind 

9 

Deg. 

AZ 

Deg. 

AZ 

mM 

6.  Z 

Ft.  MSL 

Ft 

Lb. 

12.  Z 

Ft.  MSL 

Wind 

Knots 

e 

Deg. 

Wind 

Knots 

AZ 

Deg. 

AZ 

Deg. 

A minimum  of  two  wind  points  must  be  specified.  Conditions  at  Zg  must  be  the 
first  point.  Conditions  at  Zg  must  be  the  last  point.  A maximum  of  twelve  wind 
points  may  be  specified. 
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3.3.7  SAMPLE  INPUT /OUT  PRINT 

The  following  are  copies  of  the  HP  printed  tape  for  Test  No.  3 in  Section  4. 

A.  No  Intermediate 
Altitude  Print 


PR0G#?7. 007 
3-ri  I M.  TETHER 
TEST 

3 . 000* 

RUH#1  HflX.flLT 

14000 . 000 

4000 . 000 
3.  142 
0.28  0 

25 . 000 

200. 000 

3200 . 000 

85 . 000 

WINDS 

14000 . 000 

25.0  0 0 

1 80 . 000 

1 3000 . 000 

25 . 000 

225 . 000 

1 1000 . 000 

35 . 000 

270 . 000 

113000 . 000 

40 . 000 

300. 000 

8500.000 

30. 000 

270 . 000 

4 0 0 0 . 0 0 0 
1 5 . 000 

210. 000 


ON  SURFfiCE 
ALT 

4019.508 

HT 

9980.492 

TENSION 

2951 . 693 
C. ELEV. flNG 

78.986 

LENGTH 

10080.000 

WT 

252.000 

'■,\DRRG 

38.498 

H,  I , L 

1005.997 

321 . 072 
1298.533 

fi2.T0  BLN 

39.221 

B. EL. ANG 

82.587 

S.  R 

10064.611 
CABLE  A2 

53. 198 

X +E 

821.072 

Y +N 

1005.997 

OPT. ENT 

0- NEW  PR OB 

1- NEW  “SAME  ALTS 
/WINDS 

2- LOW  CYCLES 

3- HOLD  C. LENGTH 
CHG. WINDS*HRVE 
77.007B  CARDS 
READY 

2. 000* 
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RUN# 

RUN# 

RUN# 

1 y y y 

3.  @00 

5. 000 

B. RLT/flZ 

B. RLT/flZ 

B. RLT/RZ 

13000. 000 

1 1000. 000 

8500, 000 

225.000 

270. 000 

270. 000 

3200.000 

3200. 000 

3200 . 000 

85.000 

85. 080 

85. 000 

ON  SLIRFfiCE 

ON  SURFACE 

ON  SURFACE 

RLT 

ALT 

ALT 

4011. 806 

4015. 719 

4010. 652 

HT 

HT 

HT 

8988. 1 94 

6984. 281 

4489.  34.8 

TENSION 

TENSION 

TENSION 

2976.535 

3026. 348 

3087. 994 

C.ELEV.RNG 

C.ELEV.RNG 

C.ELEV.RNG 

77.003 

"7  C'  O 7 c 
» O f O ( J 

82. 274 

LENGTH 

LENGTH 

LENGTH 

9120.000 

7080. 000 

4520. 000 

UT 

WT 

WT 

228. 000 

177. 000 

1 13. 000 

V . D R fl  U 

V. DRAG 

V. DRAG 

59.538 

46.553 

14.858 

Hj  I , L 

Hj  IjL 

H»  I j L 

1370.465 

1130.483 

519. 743 

538.72? 

75. 843 

-36. 915 

1472.550 

1133.025 

521 . 052 

flZ.TO  BLN 

HZ. TO  BLN 

RZ.TO  BLN 

66 . 460 

93.838 

85.937 

B. EL. RNG 

B. EL. RNG 

B. EL. RNG 

80 . 696 

80. 785 

83. 380 

S . R 

S.R 

S.R 

9108. 020 

7075 . 587 

4519'.  485 

CABLE  RZ 

CABLE  AZ 

CABLE  RZ 

70. 065 

8 9 . ;-!  8 9 

81 . 568 

X +E 

X +E 

X +E 

1 350 . 003 

1130. 483 

519. 743 

Y +N 

Y +N 

Y +N 

588. 128 

-75.843 

36. 915 

OPT. ENT 

0- NEW  PROB 

1- NEW  -SAME  ALTS 
/WINDS 

77. 007 


B.  Full  Print  of 
All  Intermediate 
Altitude  Data 


13501.903 
498.097 
43.578 
0. 000 
498.097 
43.578 
0.  000 


- Z,  Altitude,  Bottom  First 

- j Element 

- h 

- i 

- J 

- H 

- I 


PR0G#77. 007 
S-DIM. TETHER 
TEST 

3.000* 

RUN#1  HflX.flLT 

14000. 000 

4000. 000 
3.142 
0.280 

25.000 

500. 000 

3200.000 

85.000 

WINDS 

14000.000 

25.000 

180.000 


13000,000 

25.p00 

225.000 


500.000  - Sk, Cable  Length 
12.500  - rw.  Cable  Weight 
25.000  - Wind, knots 

202.414  - Azimuth  of  Wind 

180.000  - Sfi  + fin 

22.414  - 0!  ® 

4468.246  - R,  Reynolds  Number 
0.980  - Cd 

1.393  - q.  Dynamic  Pressure 
15.822  - E).p,  Total  Element  Drag 

1.283  - D^^r,  Vert.  Drag  Component 

14.6ii’l  — Dtt 

6.051  - Dg 


1.  181 
84. 715 
3187.593 


13004.029  - Start  of  next  element 
497.874  printout 

46. 049 
0.  949 
995.971 
89.627 
0.949 


1 1000.000 

35.000 

270.000 

10000.000 

40 . 000 

300.000 

8500.000 

30.000 

270.000 

4000.000 

15.000 

210.000 


1000,000 
25.000 
25.000 
224.819 
181. 181 
43. 638 
4532.719 
0.982 
1.418 
16.256 
1 . 084 
11.738 
11. 193 

2.099 

84.477 

3175.195 
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4095.816 
491.216 
54. 894 
75.456 
9904.  184 
976.57? 
806.570 


10000.000 
250.000 
15.319 
211 . 278 
233.964 
337.313 
3432.494 
0.980 
0.685 
7.835 
1 . 349 
7,121 
-2.977 

-0.  ;:'04 
79,058 
2955.245 

4046.725 

49.091 

5,624 

7.645 

9953.275 

982.201 

814.215 

10050. 000 
251 . 250 
15.156 
210.623 
233 . 660 
336 . 963 
3398.281 
0.980 
0.671 
0.767 
0.  134 
0.695 
-0. 296 


-0. 030 
79,040 
2954.018 


ON  SURFACE 
ALT 

4046.725 

HT 

9953.275 

TENSION 

2954.018 
C.ELEV. ANG 

79.040 

LENGTH 

10050.000 

WT 

251.250 

V.DRAG 

36,435 

H,I,L 

982.201 

314.215 
1275,799 

R2.T0  BLN 

39.658 

B. EL. ANG 

82.696 

S . R 

10034.707 
CABLE  A2 

53.630 

X +E 

814.215 

Y + N 

982.201 

OPT. ENT 

0- HEbJ  PROB 

1- NEW  -SAME  ALTS 
/bJINDS 

2- LOW  CYCLES 

3- HOLD  C. LENGTH 
CHG. WINDStHAVE 
77.007B  CARDS 
READY 
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In  the  printouts  above.  Case  A shows  the  input  data  for  a balloon  at  14,  000  ft, 
left  side  of  first  page,  followed  Immediately  with  the  surface  or  winch  data  on  the 
right  side  of  the  page.  At  the  end,  an  Option  2 or  lower  altitude  cycle  rerun  mode 
was  selected  from  4 choices  provided.  On  the  next  page.  Runs  2,  3,  and  5 are 
shown.  Run  4 is  '.'ot  shown.  Run  2,  for  example,  indicates  the  balloon  is  at 
13,000  ft  with  a wind  azimuth  of  225®.  The  same  values  of  and  6 as  in  Run  1 
were  entered  — an  approximation  since  F.^,  and  0 change  with  altitude  and  wind 
speed.  Computations  then  commence  leading  directly  to  the  printing  of  the  surface 
or  winch  conditions.  Since  an  altitude  of  8500  ft  is  the  last  altitude  above  the  sur- 
face in  the  wind  field  table  (Run  1),  the  Run  5 for  this  altitude  is  the  last  of  the 
possible  rerun  cycles.  The  program  then  terminates  with  a shorter  option  rerun 
message,  0 or  1. 

On  the  next  page.  Case  B,  the  same  balloon  problem  is  used  in  the  basic  form 
of  the  program  where  all  parameters  are  printed  at  Intermediate  altitudes  between 
the  balloon  and  the  surface.  The  altitudes  are  determined  by  the  location  of  the 
bottom -end-point  of  each  cable  element  and  are  therefore  a function  of  the  element 
length  selected  and  the  elevation  angle  of  each  element.  On  the  right  side  of  the 
page,  the  parameters  for  an  altitude  of  13,  501.  903  ft,  the  bottom  of  the  first  ele- 
ment, are  shown  followed  by  13,004.029  ft,  the  second  element.  On  the  next  page, 
the  last  two  Intermediate  altitudes  are  shown  followed  by  the  final  surface  or  winch 
printout. 

It  may  be  noted  that  while  the  same  balloon,  cable,  altitudes,  and  winds  w(3re 
specified  in  the  two  cases,  the  element  length,  K,  was  200  ft  in  Case  A and  500  ft 
in  Case  B.  As  a result,  the  intercept  surface  altitudes  differ  by  about  27  ft  for 
these  two  specific  cases.  Because  of  this  and  the  different  averaging  over  the  whole 
altitude  range,  the  other  surface  parameters  also  differ  by  small  amounts  as  will 
be  discussed  in  Section  4.  1. 
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3.4  ProCTam  No.  77.007B 


This  program,  unlike  77.007  and  77.007P,  solves  a case  where  the  cable  length 
Is  known  and  held  fixed,  while  various  wind  fields  are  introduced  and  the  balloon 
permitted  to  rise  or  fall  to  a different  equilibrium  altitude. 

3.  4.  1 SPECIAL  NOTES  FOR  OPERATION 

As  indicated  in  Table  2,  this  program  is  entered  only  after  a solution  is  found 
using  77.  007  for  a fixed  balloon  altitude,  cable,  and  wind  profile.  Entry  of  cards 
for  77. 007B  is  called  automatically  whenever  OPT.  ENT -3  is  selected  at  the  end  of 
a problem  solution  in  77.  007  (see  Section  3.  3.  2). 

With  Program  77.007B  entered,  the  storage  registers  containing  the  needed 
parameters  from  the  Program  77.007  solution  are  left  intact.  The  mathematics 
and  procedures  are  only  slight  variations  of  those  used  in  Programs  77.  007  and 
77.007P. 

A complete  solution  to  the  surface  Is  made  which  yields  a cable  length.  An 
incremental  altitude  (AZ  = 200  ft  is  built-in  but  may  be  changed  in  STEPS  0044, 

0045,  0046,  and  0047)  increase  is  then  made  followed  by  the  solution  to  a second 
cable  length.  These  cable  lengths  are  compared  to  determine  if  they  are  progressing 
in  the  proper  "direction"  towards  the  original  cable  length.  Continued  altitude  in- 
crements, either  up  or  down,  are  made  until  the  cable  length  solution  nearly  equals 
the  original  value,  A vernier  solution  with  AZ  = 20  ft  similar  to  the  K/  10  procedures 
is  then  made  for  more  precise  conditions.  This  logic  had  to  be  used  since  unlike  a 
two-dimensional  wind  field,  an  increase  in  wind  magnitude  for  example,  does  not 
always  result  in  a lower  balloon.  Wind  azimuth  changes  with  higher  wind  velocities 
on  the  cable  could  produce  a higher  balloon  condition. 

The  balloon  total  force,  P,j,  and  its  angle,  d,  are  needed  as  Inputs  into  this 
program.  However,  the  balloon  altitude  and  the  wind  at  that  altitude  are  unknown 
until  a solution  is  found.  Therefore,  an  estimate  of  Frp  and  9 must  be  made  for  the, 
as  yet,  unknown  balloon  altitude.  Some  repeat  runs  may  be  found  necessary  before 
the  balloon  altitude,  the  wind  at  that  altitude,  and  the  F.^,  and  9 values  are  all  com- 
patible, As  experience  is  gained  with  one  particular  balloon  and  cable,  these  esti- 
mations will  become  more  exact.  Repeated  runs  with  Programs  Nos.  76,003  or 
76,  005  for  a particular  balloon  at  various  altitudes  and  wind  magnitudes  can  be  used 
to  produce  a chart  of  F.p  and  9 for  any  particular  balloon  as  an  aid  in  their  selection. 
Safe  balloon  operation  would  dictate  that  in  the  original  solution  under  Program 
77,007  either;  (a)  that  the  altitude  selected  be  the  highest  allowable— envelope  full, 
balloonet  empty— and  that  any  wind  changes  permlssable  hereunder'  would  only  lower 
the  balloon,  or  (b)  that  the  proper  ballonet-fullness  effects  be  Included  in  the 
original  F,j,  and  9 used  with  the  problem  in  77,  007. 

The  wind  profile  is  entered  in  the  same  manner  as  in  Programs  77. 007  or 
77.007P  except  for  one  additional  requirement.  Because  of  the  logic  used,  one 
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positive  increment  always  moves  the  balloon  up  in  altitude  from  its  original  value 
regardless  of  whether  the  eventual  solvrtion  shows  the  balloon  higher  or  lower. 
Therefore,  the  first  wind  entry  must  be  for  an  altitude  higher  than  the  original 
balloon  altitude  by;  (a)  an  amount  sufficient  to  encompass  any  probable  final  balloon 
altitude,  or  (b)  by  at  least  200  ft  if  it  is  expected  that  the  final  balloon  altitude  will 
be  lower  than  the  original  value. 

No  prints  of  cable  parameters  at  intermediate  altitudes  between  the  balloon  and 
surface  are  provided  for  two  reasons.  Simple  print  commands  would  produce  ex- 
cessively printed  tape  since  many  interim  solutions  from  balloon  to  surface  can  be 
made  before  the  correct  solution  is  found.  Logic  allowing  printing  of  only  the  final 
solution's  intermediate  altitude  cable  parameters  requires  excessive  program  steps. 
The  balloon  altitude  and  cable  length  (or  cable  hor.  or  tension  0 if  no  solution  to  the 
surface  is  possible)  are  printed  for  each  altitude  being  tried  as  a monitoring  aid. 
When  the  correct  altitude  is  found,  the  balloon  and  surface  (winch)  conditions  are 
printed.  If  the  cable  space  position,  shape,  and  other  parameters  are  needed,  the 
resulting  balloon  altitude  and  other  input  data  can  be  used  in  Program  No.  77.007 
for  print  only  or  in  Program  No.  77.007P  for  print  and/or  plot. 

Since  incorrect  high  balloon  altitudes  are  tried  in  certain  cases  in  the  search 
before  the  correct  altitude  is  found,  the  chances  of  finding  either  zero  cable  tension 
or  horizontal  cable  conditions  above  the  ground  are  much  greater  in  operating  this 
program.  To  prevent  this  condition  from  improperly  stopping  the  program,  special 
treatment  had  to  be  evolved.  When  one  of  these  conditions  is  reached  for  the  first 
or  second  (higher)  balloon  altitude  run,  the  balloon  altitude  is  dropped  by  10  percent 
of  its  original  height  above  ground  and  the  first  and  higher  second  balloon  altitude 
calculations  are  made  again  from  this  point.  If  either  of  these  runs  yield  the  cable 
tension  zero  or  horizontal  above  ground  message,  another  10  percent  drop  is  made. 
In  this  way  the  proper  solution  is  approached  by  either  at  least  three  rising  balloon 
altitude  increments  or  any  number  of  decreasing  balloon  altitudes.  If  no  final  solu- 
tion occurs  in  this  procedure  before  the  tension  zero  or  cable  horizontal  condition 
is  signaled  on  a third  or  higher  altitude  loop,  it  can  be  properly  concluded  that  no 
solution  exists  and  such  final  message  is  given  on  the  tape.  Also  if  the  large  10 
percent  decrements  in  balloon  height  finally  place  the  balloon  at  the  surface,  the 
no-solution  message  is  given. 

3.4.2  RERUN  OPTIONS 

At  the  conclusion  of  the  printout  of  the  final  solution  of  balloon  altitude,  wind 
magnitude  and  direction  and  all  of  the  same  surface  parameters  that  are  obtained 
with  Programs  77,  007  and  77. 007 P,  a choice  of  two  rerun  options  Is  provided.  The 
option  numbers  0 and  3 have  the  same  nature  of  meaning  as  those  in  the  basic  pro- 
gram, 77,007. 
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Option  0 calls  for  a completely  new  problem  which  must  be  solved  In  77.  007. 
Hence.  If  this  option  Is  chosen  the  program  cards  for  77.007  must  be  available  for 
the  automatic  read -In  built  Into  this  program. 

Option  3 calls  for  the  cable  length  still  to  be  held  fixed  and  new  winds  to  be 
Introduced.  Use  of  this  option  then  causes  a return  to  the  start  of  77.  007B  for  re- 
entry of  F.J,,  Q,  and  winds. 


77.  007B 


3.4.3  TLOW  CHART 
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PNTtTENSion«0 


PNTj Cable 

HORlaonta 


Go  To  1077 • 


18  0m«0? 


77.  007  B 


No  _/D779\  Ytlg 
<IS  C “ 0>r^7:r 


^Is  r«17 


/^0792  \ 

Yos  \.No 

^K„7>- 


No  y'  0801  \,^^Yao 

r ClasK  >£K  ? 

V P/0838 


„ /0807  \ „ 

No  Yoa 


laiK  -kKq? 


NO  /0812\ 

— ■<Ia£K<£K„7>:- 


UOy^  Yoa 

[-<;iarK -ek^Pj^  - 


NO  //°‘"’\Yoa 


0838  I 

! 1 

Ch9  Sign  AZ 


PRKV"®b“*’® 
AZ'-  Al/10 
storo  r-1 


C'  - C + 1 

PNTi  Cablo  Longth.ZK 
Zg'  - Zg  + AZ 
Storo  Zg’  in  003,005 

OoTo  SUB/R  1480 

From  SUlVR.Clrd  Sum  Bog, 

Rnsot  K,K  Rod 

Oo  To  0258  — » — 

■>— 1 ^4 

MoioaagoiON  sURFaco  with 

AJPPROXimntoly  SANK 
LENGTH  OF  CABLE  OUT 

PNTi  BttLloon  ALTltudn  Zg 

PNTi  Balloon  WIND  Wind 

PNTi  Balloon  wind  AZimuth  Bo  "A* 

D 

PNTi  Surface  ALTitude  Z 

PNTi  HelqhT  J 

PNTi  tens ion  T 

I PNTi  Cable  ELovation  ANGlo  P 

PNTi  cable  LENGTH  1:K 

PNTi  cable  WelghT  £W 

PNTi  total  Vertical  Drag  !;d„ 

PNTi  H,I,L  H 

2 2 ^ 

L - 1-  L 

4 « arc  sin  H/L 

(UsodiTo  Polar  for  L t 4) 

AZb  -Bg  -90  - 4 

Goto  SUB/R  1687 

Prom  SUB/R,  Oriented  AZg 

PNTi  Azimuth  to  BaLloon  AZ 

B 

r - arc  tan  J/L 

PNTi  Balloon  ELevation  ANQle  e 
S.R.-L/  COB  c 

PNTi  Slant  Range  S.R 

Az^  =i:b  +Bb-1B0 

Goto  SUB/R  1687 

Prom  SUB/R,  Oriented  AZj, 

“ PNTi  CABLE  Azimuth  AZ„ 


_* 


1230 
fis  AZn 


^'11. 


Refr~450 
I 
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3.4.4  OPERATING  INSTRUCTIONS 


Four  aides  of  the  two  Program  No.  77.  007B  cards  will  have  been  loaded  after 
OPT  ENT  3 was  specified  according  to  Program  No.  77.007  instructions— Section 
3.3.4. 


KEY  STROKES  ENTRIES 

I END  1 

I FIX  I [i]>  0.  — (No.  of  decimal  places  desired) 

“contI 


(X)  (Y)  (Z) 

STOP  0.  Enter;  Test  No. 

I contI 


STOP3.  Enter;  9 

(Angle  of  (Balloon 

F.J.  to  Tot. 

Horizon,  Force  to 

deg)  Horizon, 

lbs) 

contI 


STOP  4,  Enter; 


Icont 


AZ 

(Azimuth 
of  Wind, 
deg) 


WIND 

(Wind, 

knots) 


Z 

(Altitude, 
ft  MSL) 


PRINTING 


PROG#77.  007B 
CONST.  CABLE 
LENGTH 
TEST  # 


T.  No. 


ORIG.  COND 

^B 

^S 

TT  or  Cj^ 

D 

Wt/1000 

K 

(Cable  Length)  TK^ 
EST.  NEW  FT/ANGLE 


9 

NEW  WINDS 
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Z 

WIND 

AZ 
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See  Sections  3.4.1  and  3.  4.  5 for  notes  on  first  or  highest  altitude  entry.  Stop  4*8 
will  repeat  until  last  set  of  entries  of  Zg  are  inserted.  A total  of  12  sets  of  entries 
may  be  made.  After  Zg  entry  and  CONT: 

TRIAL  ALT 


If  no  cable  “horizontal  or  tension-zero  conditions 
are  encountered,  the  progression  of  trial  altitudes 
and  their  respective  computed  cable  lengths  are 
printed  as  shown: 


SK 

Z' 

SK' 


Z" 

EK" 


When  a condition  is  met  where  the  computed  cable  length  equals  the  original  length, 
the  following  printout  is  made  including  the  abbreviated  names  and  values  of  the 
parameters  shown. 

ON  SURF.  W.  APPROX 
SAME  LENGTH  OF 
CABLE  OUT 


BLN.  ALT 

^B 

Balloon  Altitude,  ft  MSL 

B.  WIND 

Wind 

Balloon  Wind,  knots 

B.  AZ 

|3b=AZ 

Balloon  Wind  Azimuth,  Deg.  (Also  Balloon 
Pointing  Azimuth) 

S.  ALT 

Z 

Ending  Surface  Altitude,  ft  MSL 

HT 

J 

Vertical  Height,  ft 

TENS 

T 

Cable  Tension  (Winch),  lbs 

C.  EL.  ANG 

d 

Cable  Elevation  Angle  (Winch),  lbs 

LENGTH 

SK 

Cable  Length,  ft 

WT 

EW 

Cable  Weight,  lbs 

V.  D 

Ed 

w 

Total  Vertical  Drag  Comp. , lbs 

H,  I,  L 

H 

Total  H.  Distances  on  Yg  or  Y^  axis,  ft 

I 

Total  H.  Distance  on  Xg  or  X^  axis,  ft 

L 

Min,  Direct  H.  Diet,  to  balloon 

AZ  TO  BLN 

AZb 

Azimuth  Angle  to  Balloon,  deg 

B.EL.  ANG 

e 

Elevation  Angle  to  Balloon,  deg 

S.R 

SR 

Slant  Range  to  Balloon,  ft 

CABLE  AZ 

AZc 

Azimuth  Angle  of  Cable  (out  of  winch),  deg 

X 4 E 

X 

X Coordinate  to  Balloon,  ft 

y 4 N 

Y 

Y Coordinate  to  Balloon,  ft 
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If  cable -horizontal  or  tension -zero 
conditions  are  encountered,  one  of 
the  two  types  of  printoxits  shown 
will  occur  in  the  succession  of 
trial  altitudes  before  the  final 
solution  is  found  and  printed. 


C.  HOR 

TENS  0 

Z* 

Z* 

C.  HOR 

TENS  0 

Z*' 

Z*' 

SK' 

Tens  0 

z*" 

z*" 

SK" 

£K" 

Z*'" 

7*11 1 

SK'" 

EK'" 

At  this  point  the  problem  is  solved.  Two  optionsd  ways  of  continuing  are  now 
provided  by  the  following  Message  and  STOP. 

OPT  END 

0-NEW  PROB-HAVE 
77.007  CARDS 
READY 

3 -HOLD  C.  LENGTH 
CHG.  WINDS 

STOP,  Enter;  0 or  3 in  (X) 

^on3' 


If  0 is  Entered  - New  Problem— Use  this  for  a completely  new  problem  requiring 
use  of  Progrcim  No.  77. 007.  Insert  side  1 of  77. 007  Program  card  before  pressing 
the  CONT.  Continue  feeding  additional  3 sides  of  cards  and  follow  instructions  in 
Section  3.3.4.  The  permanent  storage  registers  containing  density  and  drag  coef- 
ficient constants  will  be  retained  .so  that  no  read-in  of  the  data  cards  will  be  re- 
quired-same retention  in  either  direction  bet’./een  the  two  programs. 

If  3 is  Entered  - Hold  Cable  Length,  Change  Winds— In  this  case,  the  original  cable 
length  will  be  retained,  the  summary  registers  will  be  cleared,  and  a return  to  the 
start  of  77. 007B  will  be  made.  The  first  trial  balloon  altitude  will  then  be  the  final 
balloon  altitude  solved  in  the  previous  problem. 

Note  No.  1;  In  order  to  speed  up  the  process  where  the  final  balloon  altitude  will 

be  significantly  different  from  the  starting  balloon  altitude,  an  altitude 
"guess"  value  may  be  inserted  into  Register  003  at  the  Test  Number 
STOP  for  an  original  problem  run  or  at  the  OPT  ENT  STOP  for  other 
runs. 

Note  No.  2;  If  incorrect  data  is  believed  to  have  been  entered  do  rmt  press  STOP 
END  to  restart  program.  For  correct  and  safe  clearing  of  registers, 
press  the  following; 
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77. 007B 


STOP 


LOQ.mJ 

m 

m 

0 

0 

After  OPT  ENT  message  is  printed,  enter 

1cont| 

3.  4,  b INPUT  DATA  FORM 


then: 


Test  No.  Bate:  Notes; 


ORIGINAL  VALUES 

1.  Balloon  Altitude 

= 

ft  MSL 

4.  I 

't 

= 

lbs 

2.  Wind  at  Balloon 

= 

knots 

5.  e 

deg 

3.  Azimuth  of  W,  at  Balloon  = 

degree 

6.  Cable  Length  = 

ft 

ESTIMATES 

#1 

#2 

#3 

#4 

#5 

7.  New  Balloon  Altitude 

^B 

8.  New  Wind  at  Balloon 

Wind 

9.  New  AZ  of  W at  Balloon  AZ 

1 INPUT 

1 

STOP  3 Total  Force.  F 

.j.  lbs 

Angle  of  F.J,.  0 

deg 

STOP  4 

NEW  WIND  PROFILE 

1.  Z 

6.  Z 

9.  Z 

Wind 

Wind 

Wind 

AZ 

AZ 

AZ 

2.  Z 

6.  Z 

10.  z 

Wind 

Wind 

Wind 

AZ 

AZ 

AZ 

3.  Z 

7.  Z 

11.  Z 

Wind 

Wind 

Wind 

AZ 

AZ 

AZ 

4.  Z 

8.  Z 

12.  Z 

Wind 

Wind 

Wind 

AZ 

AZ 

AZ 

Up  to  12  wind  points  may  be  specified.  First  point  must  be  a minimum  of  200  ft 

above  original  balloon  altitude  (1)  or  above  best  estimate  of  higher  balloon  altitude 

(7),  whichever  is  higher. 

Last  point  must  be  for  surface  altitude. 
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77. 007 H 

.Storuyt)  HugieitBi'H 


$ I or, AGE 

-y--"— ■«T7r--jr 

000 

00 \ l-T  T 
002  

003 

004  jS.g. 

005  J| 

006  K 

007 

008 


008 

009  

010  ~A 

on  

012  77"  <’ft  Qt> 

013  

014  Ai£.. 

0.15  oc 

016  _3:S 

017  Dtv 

018 

019  

020 

021  H 


I 

61 

J" 

Ui 

<5 

^X\l 

tc 

i. 

a 

CU<S0 

040  Vi 

041 

042 

043  W 

044 

045  "Z-i 

046 

047  /\gj 

048  Z ^ 

049  Vw- 

050 

051  g-s' 

052  V.g 

053  /4Z.; 

054  Zi, 

055  Vih 

056  AZi, 

057  g? 

058  Vt 

059  AZr 

060  Z'<? 

061  W 

062 


1 

066 

067  \/;t> 

068 

069  -J^n 

070  V// 


080  1 

081  y.2. 

002  kg  J 

083  kn  ~] 

084  ' 3 
086 

006  Rh  ] 

087  Coa..j±.. 

088  Ka  ■> 

089  Rci  I 

090  i-  5^"" 

091  _ikZ._ 

092  /6. 

093  Cog  > ^ 

094  Kg  J 

095  Rb  > 

096  Cpft  ?•  7 

097  Kr  -J 

098  ...^a.  1 

099  .SaA-.f-  ^■ 

100  /<<•/ 

101 

102  oR/(» 

103 

104  Sfitf’ZS 

105  Fr 

106  :g><  /<»8g/ 

107  C-CCONT&i 

108  © 


SP^C.TT.ooTB-*’^. Loaded  FR-oM  data  card 


77.  007B 

3.4.7  SAMPLE  INPUT /OUTPUT  PRINT 

The  following  are  copies  of  the  HP  printed  tape  for  a problem  following  Test 
No.  3 in  Section  3.  3. 7. 


ON  SURF. U. APPROX 
SAME  LENGTH  OF 
CABLE  OUT 


3.000* 

PROG#77.007B 
CONSr. CABLE 
LENGTH 
TEST  # 

3.100* 

ORIG.COND 

14000.000 

4000 . 000 
3.142 
0.280 

25.000 

500 . 000 

10050.000 
EST.NEW  FT/RNGLE 

3400.000 

79.000 

HEU  WINDS 

14500. 000 

50.000 

180.000 

12500.000 

60.000 

200. 000 

10500. 000 

65.000 

235.000 

9000.000 

55.000 

260.000 

4000.000 

35.000 

232.000 


TRIAL  ALT 

14000.000 

11500.000 

14200.000 

11750.000 

14000. 000 

11500.000 

13800. 000 

11250.000 

13600.000 

11000. 000 

13400. 000 

10700.000 

13200.000 

10450. 000 

13000. 000 

10200.000 

12800.000 

9900. 000 

12980. 000 

10150. 000 

12960. 000 

10150.000 

12940.000 

10100. 000 

12920. 000 

10100.000 

12900.000 


BLN. ALT 

12900. eet 

B.WIND 

58.006 

B.  AZ 

196.000 

S.ALT 

4026.110 

NT 

8873.890 

TENS 

3217.749 

C. EL. ANG 

51.360 

LENGTH 

10050.000 

UT 

251.250 

V.  D 

572.656 

H.  I.L 

4126.115 

1635.081 

4438.270 

A2.T0  BLN 

37.61.7 

B. EL. ANG 

63.428 

S.R 

9921.907 
CABLE  AZ 

45.084 

X +E 

2709.052 

Y +N 

3515.587 

OPT. ENT 

0-NEW  PROB-HfiVE 


77.007  CAROS 


READY 


3-HOLD  C.LEHCTH 
CHG^UINOS 


77.  007B 

The  problem  shown  above  began  with  the  entry  of  3 at  the  OPT  ENT  STOP  in 
Test  3 B,  Section  3.3.7.  A Test  No.  = 3.  1 was  assigned.  Note  that  a cable  length 
= 10,050  ft,  the  last  item  printed  under  "original  conditions",  is  the  same  as  the 
length  computed  in  the  original  problem  in  Section  3.  3.  7.  Values  of  F.p  = 3400  lbs 
and  B >=  79.  0'  were  estimated  values  for  the  new  balloon  condition. 

The  first  wind  point  introduced  here  is  at  an  altitude  of  14, 500  ft  or  500  ft  higher 
than  the  original.  Since  the  estimated  balloon  altitude  is  predicted  to  be  lower  than 
the  original,  a 200  ft  Increase  would  have  been  the  minimum  allowable.  The  new 
winds  Introduced  are  higher  in  magnitude  and  have  less  rotation  than  the  original 
case.  The  higher  magnitude  affected  the  choice  of  values  for  F'^  and  6 above. 

A series  of  trial  altitudes,  first  up,  then  progressing  downward  200  ft  at  a time, 
show  decreasing  lengths  to  an  altitude  12, 800  ft.  At  that  point,  where  the  cable  length 
is  less  than  the  original  value,  the  altitude  decrement  is  made  to  be  20  ft  but  pro- 
ceeds downward  from  the  previous  13,  OOO  ft  altitude.  Thus  12,  980,  12,  960,  etc,  are 
tried  until  at  a balloon  altitude  of  12,  900  ft,  the  equality  of  cable  length  is  established. 

The  third  column  shows  the  wind  at  the  balloon  to  be  58  knots  at  196".  The  cable 
angle  at  the  winch,  51",  is  considerably  lower  than  the  7i9"  experienced  with  the 
balloon  at  14, 000  ft.  Greater  displacements  and  other  differences  can  also  be  noted. 

Additional  procedures  might  now  be  pursued.  First  the  value  of  F.^,  and  0 could 
be  further  refined  (using  76.  003,  4,  or  5)  now  that  a "better -than -a -guess"  altitude 
is  found.  Another  runthrough  with  these  values  would  modify  the  values  of  all 
parameters  computed  in  the  first  run.  When  all  are  compatible,  some  of  the  results 
can  be  used  as  Inputs  into  77.007  for  print  only  or  77.  07P  for  print  and  plot  of  the 
complete  cable  conditions  and  geometry  between  the  balloon  and  the  surface. 
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3.5  Program  No.  77.007P 

This  version  of  the  three-dimensional  tether-cable  program  provides  a printed 
and  plotted  output.  It  requires  the  following  hardware; 

HP  Model  9810A  Calculator 
HP  No.  11210A  Math  ROM 
HP  No.  11261A  Printer-Plotter  ROM 
HP  Model  9862A  Plotter 

If  no  plot  is  desired  or  no  plotter  equipment  available,  the  use  of  Program  No. 
77.007  (Section  3.  3)  is  suggested.  Program  No.  77.007P,  however,  can  be  run 


without  the  plotter  and  with  either  a Printer-Alpha  ROM  (HP  No.  11211A)  or  the 
Printer-Plotter  ROM  by  making  the  following  substitutions  in  the  program  listing. 

Section  3,5,6; 

Step  No. 

Key 

0344 

CNT 

0345 

CNT 

0346 

CNT 

3.5,1  PLOTTING 

The  use  of  only  2 rerun  options  in  Program  No,  77, 007 P provides  more  pro- 
gram stops  than  Program  No.  77.007.  This  permits  the  plotting  of  the  following 
six  curves  with  symbols  as  shown  on  one  piece  of  ~ 10x16  in,  paper. 


(a) 

Symbol 

+ 

Parameters 

Z-H 

(d) 

Symbol 

T 

Parameters 

Z-T 

(b) 

« 

Z-I 

(e) 

0 

z-e 

(c) 

, 

H-I 

(f) 

X 

Z-q 

Before  computation  and  plotting  of  data  commence,  various  events  take  place 
to  prepare  the  plot  paper.  Axes  for  curves  a,  b,  and  c are  drawn  as  follows; 
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In  addition  a fifth  line  is  drawn  along  the  top  edge  of  the  paper,  3 -in.  above  the 
upper  I -axis  for  the  tension,  d.  During  the  drawing  of  axes,  tick  marks  are  added 
at  user  specified  intervals.  No  tick-marks  for  the  angle  6 are  provided;  the  scale 
for  6 is  always  10®  per  inch  or  90°  is  9-in.  to  the  right  of  the  Z-axis.  No  tick- 
marks  for  q are  provided;  the  scale  for  q is  always  2-lbs/ft  per  inch. 

Three  constants  are  entered  at  STOP  5 in  order  to  scale  the  axes  and  set  tick- 
mark  spacing.  The  three  constants  are; 

N.  The  interval  in  full-scale  feet  between  tick-marks  on  the 
Z,  H,  and  I axes.  Zg  OR  SURFACE  ALTITUDE  is  located 
at  the  origin  of  the  Z-H  or  Z-I  set  of  axes  with  Zg  located 
at  the  top  of  the  Z-axis.  Since  the  available  length  of  the 
Z-axis  is  10-in. , the  magnitude  of  the  altitude  excursion, 

Zg-Zg,  will  determine  the  scale  of  the  Z,  H,  and  I axes 
and  the  tick -mark  intervals.  Therefore  in  the  case  of  a 
high  altitude  balloon  with  light  winds  and  a small  amount  of 
cable  displacement,  the  H - I plot  may  be  quite  small. 

M.  The  Interval  in  pounds  between  tick-marks  on  the  tension 
axis. 

P.  A factor  used  for  scaling  the  tension  axis.  The  value  of 

p 

the  balloon  total  force,  P.^,,  rounded  to  the  nearest  10 
pounds  is  located  6 -in.  to  the  right  of  the  Z axis.  Tick- 
marks  are  drawn  at  M intervals  with  one  at  the  Z axis  where 
TsO.  The  rounding  of  F.j,  above  is  done  only  for  establishing 
the  scale  of  the  T axis— the  computations  are  unaffected.  If 
for  example,  F = 7075,  P =>  2,  and  M = 1000  the  process 
would  round  F.^,  to  7 100  lb  and  divide  by  6 to  get  the  scale 
factor.  Using  the  scale  factor,  the  1000-lb  intervals 
would  be  converted  to  the  proper  linear  dimension  for 
marking  on  the  paper.  When  plotted,  the  exact  F.p  point 
would  appear  just  to  the  right  of  the  7000-lb  tick -mark. 

The  plotting  paper  will  then  appear  as  shown  on  the  following  page  (without 
notations)  before  any  computations  or  plotting  commence;  the  upper  axis  to  the  left 
of  the  Z-axis  represents  imaginary  negative  values  of  T and  6 ^nd  can  be  ignored. 
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The  spacing  of  the  axes  favors  conditions  of  positive  H and  I values  although 
negative  values  of  lesser  magnitude  can  be  accommodated.  The  H-I  plot  is  really 
a Xg  - Yg  plot  looking  down  on  the  balloon  and  cable  from  above.  The  balloon, 
located  at  the  origin  of  the  H and  I axes,  is  pointing  down  the  paper  in  the  positive 
H direction  the  same  as  in  Figure  4.  Similarly,  the  wind  at  the  balloon  is  pointed 
up  the  paper.  Knowing  this  balloon  wind  azimuth,  it  is  then  simple  to  draw  a North 
reference  line  on  the  plot  with  a protractor  if  desired.  Similarly  the  balloon  in  the 
vertical  plane  is  located  at  the  top  of  the  Z-axls  at  Zg  and  is  pointing  to  the  right. 

If  a two  dimensional  wind  problem  were  plotted,  the  Z-I  curve  would  be  hidden  in 
the  Z-axls  (1=0)  and  the  H-I  curve  would  be  hidden  in  the  vertical  H-axis  (1=0). 

The  number  of  points  plotted  for  each  of  the  six  cuiwes  will  equal  the  number 
of  cable  elements  required  to  reach  the  surface.  If  any  plotted  point  should  be 
called  that  exceeds  the  10  X 15 -in.  boundaries  of  the  paper,  it  will  not  be  plotted 
but  the  other  points  will  still  be  plotted.  See  also  Section  3.  5.  7. 

3.5.2  RERUN  OPTIONS 

When  a problem  is  solved,  the  cable  end  conditions  printed,  and  the  plotter  is 
stopped,  the  program  is  designed  to  offer  the  user  a choice  of  methods  with  which 
to  proceed  to  the  next  problem.  A message  is  printed: 

OPT.  ENT  (OPTION,  ENTER) 

followed  by  a choice  of  two  numbers  and  a STOP, 

When  a completely  new  problem  requiring  reentry  of  different  altitudes,  winds, 
and  cable  parameters  is  next  to  be  run,  the  number  0 is  entered.  The  program  will 
clear  all  but  the  storage  registers  loaded  from  the  data  cards  and  cycle  back  to  the 
initial  printing  of  the  program  number  and  title. 

The  normal  runs  (printed  RUN(/1MAXALT)  provide  a cable  solution  (a)  for 
the  balloon  at  a specific  altitude  subjected  to  the  winds  specified  in  the  first  group 
of  wind  entries  and  (b)  for  the  cable  subjected  to  a variety  of  wind  conditions 
specified  at  the  lower  altitudes  down  to  the  surface.  One  practical  and  sometimes 
limiting  problem  is  concerned  with  raising  or  lowering  the  balloon  through  the  same 
specified  wind  field.  Rather  than  a reentry  of  all  variables  required  at  each  of 
many  lower  altitudes,  a method  is  provided  to  simplify  a series  of  runs  at  decreas- 
ing altitudes. 

If  the  number  2 is  entered  at  the  above  STOP,  the  program  will  clear  the  sum- 
mation registers  and  retain  the  surface  altitude,  cable  specifications,  and  wind 
field,  and  setup  RUN  No.  2 with  the  balloon  at  the  second  lower  altitude  in  the  wind 
field  table.  Because  the  balloon  total  force  and  angle  at  this  new  altitude  will  differ 
from  the  previous  balloon  altitude,  only  one  STOP  is  needed  for  those  two  entries. 

A full  solution  is  then  provided  for  this  condition.  At  its  end.  the  program  loops 
to  RUN  No.  3 and  places  the  balloon  at  the  third  lower  altitude  in  the  wind  field 
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table.  These  automatic  loops  continue  until  the  surface  is  reached.  At  that  point 
no  OPT.  END  message  is  printed.  The  program  cycles  back  to  a restart  as  if 
Option  O were  selected. 

When  considering  the  use  of  this  lowering  altitude  cycle  option,  it  is  necessary 
to  initially  enter  a wind  field  table  which  specifies  a sufficient  number  of  altitudes 
to  provide  a complete  performance  analysis.  Altitudes  having  particularly  strong 
winds,  of  concern  in  balloon  performance,  would  of  course  be  included  in  the  initial 
definition  of  the  wind  field  acting  on  the  cable.  Similarly,  it  is  necessary  to  have 
the  values  of  the  balloon  total  force,  P.j„  and  its  angle,  0,  available  for  entry  at 
each  altitude  as  indicated  in  the  INPUT  DATA  FORM.  The  use  of  Program  No. 

76.  003  or  No.  76.  005,  Reference  1.  for  each  of  the  altitudes  in  the  wind  field  table 
will  provide  the  values  of  F.j,  and  0. 

When  using  the  plotter  with  the  OPT  2 mode,  it  should  be  noted  that  no  new 
axes  or  tick  marks  will  be  drawn  for'each  of  the  balloon  altitudes.  If  the  original 
graph  from  the  maximum  balloon  altitude  solution  is  left  on  the  plotter  and  only  2 
or  3 lower  altitude  solutions  were  involved,  pen  color  changes  might  help  to  keep 
each  set  identifiable.  All  scales  except  the  tension  scale  remain  constant  for  the 
lower  altitude  runs.  The  tension  scale  may  or  may  not  remain  fixed  depending  on 
the  value  of  the  newly  entered  at  each  balloon  altitude.  The  original  values  of 
M and  P entered  at  STOP  5 are  retained  and  operate  on  the  new  F.p.  If  F rounds 
to  the  same  number  as  In  the  RUN  No.  1 for  the  maximum  altitude,  then  the  scale 
remains  the  same  and  the  Z-T  curve  may  be  referenced  to  the  tension  axis  and 
tick  marks. 
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3.  5.  4 OPERATING  INSTRUCTIONS  AND  NOTES 
KEY  STROKES  ENTRIES 


1 RUN 

1 END 

FIX 

[El  0 - 

^ fil  — (Number  of  decimal  places  desired) 


LOAD  I 


Insert  Side  1 of  Program  Cards 


Continue  inserting  Sides  2,  3,  and  4 of  Program  Cards 
Insert  Side  1 of  Data  Cards 


I END  I 
CONT 


(Following  heading  is  now  printed) 
PROG  # 77.007P 
3 -DIM.  TETHER 
TEST 


(X)  (Y)  (Z) 

At  STOP  0-0-0,  Enter;  Test  No/I.  D.  

I CONT 


Following  printing  of  the  test  number  or  I.  D. , Side  1 of 
the  Data  Card  will  load. 


Continue  inserting  Sides  2 and  3 of  Data  Cards 


Note  that  loading  of  Program  and  Data  Cards  will  only  be 
required  one  time  until  machine  is  turned  off. 


(Following  is  now  printed) 


RUN  #1 

At  STOP  1-1-1,  Enter; 


ICONT  I 

At  STOP  2-2-2,  Enter; 


fcONT  1 

At  STOP  3-3-3,  Enter: 
|CONf~l 


17  or  Zg  (ft,  MSL)  Zg  (ft.  MSL) 

Enter  n to  use  built-in  cylinder  C^.  variation 
or  enter  Cq  value  which  will  be  held  constant 
throughout  program  run. 

OorK*  Wt/ 1000  ft  Diam.  (in. ) 

(Element  (lbs)  (Cable) 

Length)  (Cable) 

* - Z„ 

Enter  0 to  set  K = — or  enter  K in  ft. 


$ (deg) 
(Angle  of  F„ 
to  Horizonf 


F (lb) 
(Balloon  Tot. 
Force) 


WINDS  (Printed) 
At  STOP  4-4-4,  Enter; 


Azimuth  of  Wind  (knots)  Z_  (ft,  MSL) 

Wind  (deg) 


This  entry  must  be  for  conditions  at  balloon  altitude. 
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(X)  (Y)  (Z) 

|C0NT~ 

At  STOP  4-4-4,  Enter:  Azimuth  Wind  Z 

Stop  4's  will  repeat  until  the  last  set  of  entries  for  the  surface 
condition  at  Zg  are  inserted.  A total  of  12  sets  of  wind  entries 
may  be  made  including  the  balloon  and  surface  conditions. 

CONT 


The  above  sets  of  entries  will  be  printed  out  in  groups  when  CONT  is  struck  after 
STOPS  1.  2,  3,  and  4.  The  entered  or  computed  value  will  be  printed  in  the  case 
of  K. 


(If  no  plot  is  wanted,  see  Section  3.  5) 

The  plotter  routine  next  becomes  operative  requiring  the  following; 
Install  plot  paper  with  limits  set  to  10  X 15 -in. 

At  STOP  5-5-5,  Enter:  P M N 


CONT 


(Rounding  Factor 
lOP-lb,  for 
Tension  Scale) 


(Tension  Tick 
Mark 

Interval,  lb) 


(Spacial  Tick 
Mark 

Interval,  ft) 


Axes  and  tick-marks  will  now  be  drawn  by  the  plotter  per  Section  3.5.  1. 
Note  the  scale  of  6 is  always  IC/in.  with  0 = 0 at  the  Z-axis  and  the  scale 
of  q is  always  2 Ib/ft^  per  in. , q = 0 at  the  Z-axis. 


The  program  now  begins  computation  starting  at  the  balloon  and  works  downward 
one  element  at  a time  to  the  surface.  Twenty-three  parameters  are  first  printed 
then  the  following  six  points  are  plotted,  all  representing  the  conditions  at  the 
bottom  end-point  of  each  element.  The  six  points  are  plotted  with  the  following 
symbols; 

+ Z-H 

Z-I 

0 Z-0 

T Z-T 

X Z-q 

H-I 

Nols:  To  avoid  the  printout  of  any  or  all  of  the  twenty -three  parameters  shown 
below,  replace  PNT  with  CNT  at  the  associated  program  step  numbers. 


Step  No. 

0395 

Z 

Altitude,  ft  MSL 

0397 

j 

Vert.  Distance,  top  to  bottom  of  element,  ft 

0420 

h 

Horiz.  Distance,  parallel  to  balloon  axis,  ft 

0431 

i 

Horiz.  Distance,  perpendicular  to  bain,  axis,  ft 

0439 

J 

Total  Vert.  Distance,  balloon  to  bottom  of  element,  ft 

0443 

H 

Total  Horiz.  Dist.  , sh,  balloon  to  bottom  of  element. 

0447/8 

I 

Total  Horiz.  Dist. , Zi.  balloon  to  bottom  of  element, 

129 


77.  007P 

Step  No. 

0452 

ek 

Total  Element  (Cable)  Length,  ft 

0462 

IW 

Total  Element  (Cable)  Weight,  lb 

0563 

Wd 

Wind  Velocity  at  bottom  of  element,  knots 

0665 

AZ 

Wind  Azimuth  at  bottom  of  element,  deg 

0667 

Azimuth  of  Vert.  Plane  Containing  the  Element,  deg 

0688 

a 

Wind  Incidence  Angle  on  the  Element,  deg 

0812 

R 

Reynolds  Number 

0994 

Drag  Coefficient 

1011 

q 

'J 

Dynamic  Pressure,  Ib/ft  , Based  on  Velocity  Normal  to 
Element 

1021 

Dt 

Total  Element  Drag,  lb 

1060 

Vertical  Drag  Component,  lb 

1072 

°H 

Horiz.  Drag  Comp,  in  Vertical  Plane  of  the  Element,  lb 

1086/7 

DS 

Horiz.  Drag  Comp,  to  Vertical  Plane  of  the  Element,  lb 

1114 

^1 

Horiz.  Rotation  of  Tension  Vector  at  bottom  end  or 

Horiz.  Rotation  of  the  nexrt  element's  vertical 
plane,  deg 

1148 

1206 

^1 

Pitch  Angle  Downward  of  Tension  Vector  at  bottom  end  or 
Elevation  Angle  of  next  element  above  the  horizon,  deg 

1161/2 

1216/7 

^1 

Tension  at  bottom  end  or  at  top  end  of  next  element,  lb 

Groups  of  the  values  of  the  above  23  parameters  will  continue  to  be  printed  for  points 
down  the  cable  until  one  of  the  following  conditions  is  encountered; 

(a)  The  cable  reaches  the  earth's  surface  at  7g— the  winch  location, 

(bl  The  tension  becomes  rero, 

(c)  The  cable  becomes  horizontal. 

(In  (a),  the  computational  techniques  used  do  not  yield  a precise  Zg  condition.  The 
final  Z will  be  higher  than  Zg  by  an  amo\int  less  than  K sin  0/ 10,  usually  no  more 
than  a few  feet.  ) 

The  final  printout  includes  the  abbreviated  names  and  values  of  the  following  param- 
eters. They  describe  the  conditions  at  the  winch  if  condition  (a)  is  attained  or  at 
the  cable  lower  end  which  is  above  the  surface  if  conditions  (b)  or  (c)  are  indicated. 


(a)  ON  SURFACE 


ALT  Z 

HT  J 

TENS  T 

C.ELEV.ANG  0 
LENGTH  IK 

WT 


or  (b)  TENSION.  0 or  (c)  CABLE  HOR 

Altitude,  ft 
Vertical  Height,  ft 
Cable  Tension,  lb 

Elevation  Angle  of  Cable  above  Horizon,  deg 
Cable  Length,  ft 
Cable  Weight,  lb 

Total  Vertical  Drag  Component,  lb 
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H.  I.  L 


AZ.  TO  BLN 

B.  EL.  ANG 
S.R 

C.  AZ 
X + E 
Y + N 


II 

I 

L 


AZ 


SR 


Tot.  Iloriz.  Distance  along  Yg  or  Y^  axis,  ft 
Tot,  Iloriz.  Distance  along  Xj,  or  X^  axis,  ft 
Min.  Direct  Horizontal  Dist.  to  Balloon,  ft 
Azimuth  Angle  to  Balloon,  deg 
Elevation  Angle  to  Balloon,  deg 
Slant  Range  to  Balloon,  ft 
Azimuth  Angle  of  Cable  {Out  of  Winch),  deg 
X Coordinate  to  Balloon,  ft 
Y Coordinate  to  Balloon,  ft 
At  this  point  the  initial  problem  entered  is  solved  with  printing  and  plot  completed. 
If  this  was  an  initial  run  (RUN  #1  MAX  ALT),  the  following  is  printed  and  STOP 
provided  to  permit  2 optional  ways  to  rerun  the  program. 


AZ 


C 

X 


OPT.  ENT 
0 -NEW  PROB 
2 -LOW  CYCLES 


At  STOP  , Enter;  0 or  2 in  (X) 

IcontI 


If  0 is  Entered  - New  Problem— Use  for  completely  new  problem.  All  but  the 
permanent  storage  registers  containing  density  and  drag  coefficient  constants  will 
be  cleared,  program  returns  to  start  with  reprint  of  Number  and  Title. 

If  2 is  Entered  - Lower  Altitude  Cycles— Use  when  analysis  is  desired  with  the 
balloon  lowered  to  each  of  the  altitudes  specified  in  the  wind  field  table.  Summary 
storage  registers  are  cleared  and  the  program  returns  to  start  with  reprint  of 
Number  and  Title.  Instead  of  RUN#1  the  following  is  printed: 

RUN  # 


B. ALT/AZ 


2 


Program  then  goes  to  STOP  3 for  entry  of  balloon  total  force  and  angle  at  this  new 
altitude  Zg  which  is  the  second  altitude  originally  entered  in  the  wind  table.  The 
program  will  then  make  a complete  solution  to  the  surface  (or  to  zero  tension  or 
horizontal  cable)  for  this  new  balloon  altitude.  However  at  the  end  of  RUN  #2,  no 
option  is  provided  since  the  program  cycles  on  to  the  third  altitude  point  in  the  wind 
field  table,  sets  up  RUN  i'3  and  goes  to  STOP  3 again  for  entry  of  the  two  balloon 
parameters  at  this  new  altitude.  Thus  a solution  is  provided  for  each  altitude  in 
the  wind  field  table  except  for  Zg.  When  Zg  is  detected,  the  program  will  terminate 
the  same  as  if  0 were  entered  at  OPT  ENT  and  return  to  the  start  with  a reprint  of 
Number  and  Title.  Refer  to  Section  3.5.  1,  for  plotting  notes  with  this  OPT  2. 
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Note  No.  I - If  Inooi'k'oot  dutu  la  believed  to  have  been  oiUorod,  do  not  proas 
STOP  KNDto  roatHii*t  program.  For  eorroet  and  safe  clearing  of 
reglatera  proas  foUowlngi 
STOP 
QO  TO 
1 
S 
0 
0 

After  OPT  ENT  message  Is  printed,  enter  0,  then: 

CONT 

Note  No.  2 - Proof  of  solutions  can  be  made  by  summation  of  the  vertical  forces. 

Summary  of  the  horizontal  forces  is  not  possible  because  summation 
along  two  fixed  horizontal  axes  could  not  be  handled  within  the 
available  space. 

Therefore;  F.p  sin  6 « T^  sin  0 ^ +ED^ 
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3.  5.  5 INPUT  DATA  FORM 

Tost  No, ! Date:, Notes: 


1 input  77.007P  1 

STOP  NO, 

ITEM 

VALUE 

1-1-1 

MAX.  BALLOON  ALTITUDE  Zg 

SURFACE  ALTITUDE  Zg 

CABLE  Cg  or  v for  Internal  Comp. 

Ft.  MSL 

Ft.  MSL 

2-2-2 

CABLE  DIAMETER  D 

CABLE  WEIGHT  per  1000  ft 

CABLE  ELEMENT  LENGTH,  KORO 
for  K = (Zg  - Zg)/100 

Inches 

Lb. 

Ft, 

3-3-3 

BALLOON  TOTAL  FORCE  F,j, 

ANGLE  OF  TOTAL  FORCE  6 

Lb. 

Deg. 

4-4-4 

WIND  PROFILE 

For  Stop  3 
Opt.  Lower 
Alt.  Cycles 

For  Stop  3 

Opt.  Lower 

Alt.  Cycles 

1.  Z„  Ft.  MSL 

is 

wind  Knots 

AZ  Deg. 

7,  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j,  Lb. 

9 Deg. 

2.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j.  Lb. 

9 Deg. 

8.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

Lb. 

9 Deg. 

3.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j,  Lb. 

6 Deg. 

9.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F^  Lb. 

9 Deg. 

4.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg 

F,j.  Lb. 

9 Deg, 

10.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j,  Lb. 

9 Deg. 

5.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j,  Lb. 

9 Deg. 

11.  Z Ft,  MSL 

Wind  Knots 

AZ  Deg. 

F^  Lb. 

9 Deg. 

6.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

F,j.  Lb. 

9 Deg. 

12.  Z Ft.  MSL 

Wind  Knots 

AZ  Deg. 

A minimui 

first  point 
points  maj 

n of  two  wind  points  must  be  specified.  Conditions  at  Zg  must  be  the 
Conditions  at  Zg  must  be  the  last  point.  A maximum  of  twelve  wind 
be  specified. 

5-5-5 

PLOT  CONSTANTS 

Special  Tick-Mark  Intervals,  N 
Tension  Tick-Mark  Intervals,  M 
Tension  Rounding  Factor,  P 

Ft. 

Ft. 
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V7.  007P 


storage  Registers 


040  Vi 

041  A^i 

042  ga. 

043  W 

044  Aix. 

045  'Z.1 
04G  Vi 

047  Ajj 

048 

049  VW- 

050 

051 

052  Vs 

053  >4Zs 

054  Zc 

055  Vo 

056  Aiii 

057  g? 

058  W 

059  AEr 
OGO 

061  V? 

062 


065 

066  Z/c 
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3.  5.  7 SAMPLE  INPUT/OUTPUT  PRINT  AND  PLOT 

The  following  are  copies  of  the  HP  Printed  Tape  for  Test  No.  5 in  Section  4. 

A.  No  Intermediate 
Altitude  Print 


PROG#??. 007P 
3-IiIM.  TETHER 
TEST 

5.000* 

RUN#1 

10000.000 

4000.000 
3.  142 

1.250 

20.000 

180.000 

3200.000 

85.000 

WINDS 

10000. 000 

25.000 

180.000 

8500.000 

40.000 

225.000 

7000. 000 

50.000 

270.000 

6000.000 

68.000 

300.000 

4000.000 

20.000 

315.000 


ON  SURF 


ALT 

4003.890 

HT 

5996.110 

TENS 

3271 . 459 

C.ELEV.flNG 

35.932 

LENGTH 

7290.000 

UT 

145.800 

V.  D 

1 122.263 

Hj  I)  L 

1476.213 

3000.984 

3344.415 

flZ.TO  BLN 


B.  EL.  FING 
S.  R 

C.  flZ 
X.  +E 
Y +H 


63.807 
60.849 
6865.745 
85.548 
3000.984 
1476.213 


OPT. ENT 
0-NEW  PROB 
2-LOU  CYCLES 


156 


RUN# 

£.£itiO 

B. RLT/HZ 

SSOe. 000 

225.000 


3200.000 

5.000 


CRBLE 

ALT 

HOR 

HT 

8089.653 

TENS 

410.34? 

3209.964 
C. ELEV. RNG 

LENGTH 

-0. 136 

UT 

8280.000 

V.D 

165.600 

H»  IsL 

120. 934 

flZ.  TO 

326? . 066' 
2.584 
826?. 066 
BLN 

45.018 

B.  EL. flNG 

S.  R 

2.842 

C.flZ 

8 2 f’  ? . 244 

X +E 

45. 0?4 

Y +N 

584?. 525 

5843. 8?1 

RUN# 

RUN# 

3.000 

4,000 

B.RLT/RZ 

e.RLT/R2 

?000.O0O 

6000,000 

270.000 

300.000 

3200.000 

3200.000 

85.000 

85.000 

ON  SURF 

ON  SURF 

ALT 

RLT 

4007.570 

4005.953 

HT 

HT 

2992.430 

1994.047 

TENS 

TENS 

3171.082 

3174.239 

C. ELEV. RNG 

C. ELEV. RNG 

56.523 

68.224 

LENGTH 

LENGTH 

32.58. 000 

2070;,  000 

UT 

UT 

65.160 

41 . 400 

V.D 

V.D 

477.637 

198.685 

H»  I,L 

H»  I,L 

1 113,406 

520.478 

369.571 

28.087 

1 173. 139 

521 . 235 

flZ.TO  BLN 

flZ.TO  BLN 

108.362 

123.089 

B. EL. RNG 

B. EL. RNG 

68.593 

75.351 

S.R 

S.R 

3214. 170 

2061 . 046 

C.R2 

C.  RZ 

113.851 

124. 851 

X +E 

X +E 

1 113.406 

436. 703 

Y +N 

Y +N 

-369. 5?1 

-284.563 
PR0G#77. 007P 
3-DIM. TETHER 

TEST 

77,  007P 


l! 


4 13  m ( 4 ti,  Cl 
U3,5?^ 
M45 
I 4 . S I S 
§'3(38.  ?mK3 
147’8.CI?§ 
i'38^.,  4«3C< 

?^?£.Cic3Ci 

148.440 

314.0'?2 

3^>S.4'3i: 

4'3. 

ir84^.i:c«r 
I . so  Cl 
0. 8<»8 
S . Cl  1 '3 
1 . 088 
1 . n? 

1.804 


0. 058 
38. 

32;’ 1 . 8'i>8 
4003. 830 

10. 8n 

1 . I 38 
14.528 

5338. I 10 
14?8.213 
3000.384 

■’230.000 
145.800 
20.078 
314.371 
285.520 
43.450 
17375.023 
1.200 
0,880 
1.373 
1 . 041 
1 . 034 
1 . 273 

0.028 

35.332 

3271.453 


OH 


ftLT 

HT 

4003,830 

Tens 

5338. lie 

3271.453 

C.ELliV.ftNG 

LENGTH 

38.332 

iiir 

7230.000 

V.D 

145.800 

H.  I .L 

1 122.283 

1478.213 

3000.384 

3344.415 

82.  TO 

BLN 

83,807 

e, EL. fiHG 

8.R 

80.843 

C.flZ 

8885.745 

.S 

85.548 

Y +N 

3000.384 

1478.213 

OPT. ENT 

0-NEW 

PROB 

2-LOlJ 

CYCLES 
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In  the  printouts  above.  Case  A shows  the  input  data  for  a balloon  at  10,  000  ft, 
left  side  of  the  first  page  which  is  followed  immediately  with  the  surface  or  winch 
data  on  the  right  side  of  the  page.  At  the  end,  an  Option  2 or  iower  altitude  cycle 
rerun  mode  was  selected  from  the  two  choices  provided.  On  the  next  page,  Runs 
2,  3,  and  4 are  shown.  Run  2,  for  example,  indicates  the  balloon  is  at  8500  ft  with 
a wind  azimuth  of  225“.  The  same  values  of  F.p  and  6 as  in  Run  1 were  entered; 
an  approximation  since  F.^,  and  9 change  with  altitude  and  wind  speed.  Computation 
then  commences  leading  directly  to  the  surface  or  winch  conditions.  Run  3 auto- 
matically follows.  Since  6000  ft  is  the  last  altitude  above  the  surface  in  the  wind 
field  table  (in  Run  1),  the  Run  4 for  this  altitude  is  the  last  of  the  possible  rerun 
cycles.  The  program  therefore  terminates  by  going  back  to  the  start,  reprints  the 
title,  and  is  ready  for  a new  problem. 

On  the  next  page.  Case  B,  the  same  balloon  problem  is  used  in  the  basic  form 
of  the  program  where  all  parameters  are  printed  at  intermediate  altitudes  between 
the  balloon  and  surface.  The  altitudes  are  determined  by  the  location  of  the  bottom 
end-point  of  each  cable  element  and  are  therefore  a function  of  the  element  length 
selected  and  the  elevation  angle  of  each  element.  On  the  right  side  of  the  page,  the 
parameters  for  an  altitude  of  9820. 685  ft,  the  bottom  of  the  first  element,  are  shown 
followed  by  9641.  584  ft  for  the  second  element.  On  the  next  page,  the  last  two 
intermediate  points  are  shown  before  the  final  surface  or  winch  condition  printout. 
These  final  winch  figures  agree  with  those  in  Case  A,  Run  1,  since  all  input  data 
for  both  cases  was  identical. 

The  plot  for  this  problem  is  then  shown.  All  notations  are  added;  none  are 
produced  by  the  program  except  the  point  symbols  noted  in  Section  3.  5,  1.  Tick- 
mark  intervals  of  500  ft  for  all  spacial  dimensions  and  500  lb  for  the  tension  were 
selected  as  part  of  the  input  data,  A tension  rounding  factor  of  3 made  the  initial 
F.J,  round  from  3200  lb  to  3000  lb  thereby  placing  the  3 000- lb  tick  mark  in  the  middle 
of  the  upper  (tension  axis).  The  3200  lb  starting  value  may  be  seen  to  lie  at  the 
correct  location  to  the  right  of  the  3000-lb  tick. 

Because  of  the  unusual  cable  properties  selected  to  illustrate  large  cable  rota- 
tion (cable  diam  = 1.  25  in,  and  weighing  only  20  lb  per  1000  ft),  the  tension  tends 
to  remain  constant  and  then  increases  above  its  starting  value  as  it  approaches  the 
winch.  In  the  more  typical  heavy  cable  used  in  tethered  balloon  work,  the  tension 
usually  decreases  with  decreasing  altitude. 

The  winds  changing  from  south  past  westerly  moving  down  from  the  balloon 
produced  a large  amount  of  cable  turning  as  may  be  noted  in  the  H-I  plot  which  is 
really  an  Xg-  plot  looking  vertically  down  from  above  the  balloon.  Due  to  the 
input  value  of  wind  azimuth  at  the  balloon,  180°,  the  north  direction  is  vertical 
along  the  H-axis  toward  the  top  of  the  paper.  The  balloon,  always  pointing  down 
the  paper  at  the  intersection  of  the  H-I  axes,  is  "looking"  into  a south  wind.  The 
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77.  007P 

westerly  winds  on  the  cable  below  the  balloon,  coming  from  the  left,  "turn"  the 

cable  in  that  direction.  The  end  of  the  H-I  plot  represents  the  winch  location.  Final 

values  of  location,  azimuth  of  the  balloon  from  the  winch,  and  cable  out  azimuth 

may  be  scaled  from  the  plot  and  seen  to  agree  with  the  printed  figures.  The  two 

vertical  projections  of  the  cable  position,  Z-H  and  Z-1,  can  be  seen  to  be  consistent 

with  the  H-I  plot  since  all  three  are  at  the  same  scale. 

The  cable  elevation  angle,  6 , decreases  from  85°  at  the  balloon  to  35.  9°  at  the 

winch.  In  the  full-size  10-in.  X 15-in.  plot,  6 = 90°  is  9 in.  to  the  right  of  the 

Z a.'cis.  Knowing  that  the  two  vertical  axes  are  8 in.  apart  on  the  full  scale  plot,  it 

is  possible  to  ascertain  the  scale  of  d on  any  plot  reproduced  to  a smaller  size. 

2 

Similarly,  since  the  scale  of  q is  2 Ib/ft  per  in.  on  a 10  in.  X 15  in.  plot,  the 
scale  on  a plot  of  reduced  size  may  be  ascertained  by  reference  to  the  distance  be- 
tween the  two  vertical  axes. 

Note  that  the  dynamic  pressure,  q,  is  based  on  the  wind  velocity  normal  to 
each  cable  element  and  is  therefore  a function  of  the  wind  velocity  squared,  the 
atmospheric  density,  the  elevation  angle,  6 , and  the  wind  incidence  angle,  a . In 
this  particular  example,  q reaches  a maximum  at  6000  ft  where  the  wind  is  a 
maximum.  If  for  example,  the  cable  elevation  angle  at  this  point  were  20°  instead 
of  an  apparent  52°,  the  dynamic  pressure  would  be  reduced  due  to  the  smaller  wind 
vector  normal  to  the  cable. 


4.  PROBLEM  SOLUTIONS 

Some  examples  of  the  solution  to  various  tethered  balloon  cable  problems  will 
serve  to  indicate;  (1)  the  ability  to  handle  some  of  the  intricacies  involved  in  working 
with  azimuth  angles  in  a three-dimensional  solution,  and  (2)  the  usefulness  in  ad- 
vanced design  efforts.  Discussions  here  on  the  use  of  Program  No.  77.  007  apply 
equally  to  No.  77.  007P. 

4.1  Problems  for  Testing  Program  77.007  and  77.007P  Operation 

4.  1.  1 180°  AZIMUTH  AMBIGUITY -TESTS  1 AND  2 

A.  Test  1 

A sample  problem  used  in  Reference  1,  Program  No.  76.  006  will  be  utilized 
here  as  Test  No.  1.  In  this  example  of  a two-dimensional  problem,  the  wind  was 
reversed  180°  in  direction  between  altitudes  of  10,  000  and  8000  ft  by  changing  the 
sign  of  the  wind  magnitude  in  the  wind  profile  entries  in  76.  006,  In  the  case  of  this 
type  of  two-dimensional  problem  operated  in  the  Program  No,  77.  007,  the  same 
reversal  was  first  entered  by  a change  in  azimuth  as  shown  before. 
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User  Entries 

76. 006 

77. 007 

Balloon  Altitude,  ft  MSL 

14,  000 

14, 000 

Surface  Altitude,  ft  MSL 

4,  000 

4,000 

Internal  Cj^  Computations 

n (Yes) 

TT  (Yes) 

Cable  Diameter,  in. 

0.  28 

0.  28 

Cable  Weight/ 100  ft,  lb 

25.  0 

25.  0 

Element  Length  (K),  ft 

500,  0 

500.  0 

Balloon  Total  Force,  lb 

1385. 0 

1385. 0 

Total  Force  Elev,  Angle,  deg 

79.4 

79.  4 

Wind  Field 

Z 

14,  000 

14,000 

Wind,  knots 

25 

25 

Azimuth,  deg 

180 

Z 

10,  000 

10,  000 

Wind,  knots 

60 

60 

Azimuth,  deg 

— 

180 

Z 

8,  000 

8,000 

Wind,  knots 

-15 

15 

Azimuth,  deg 

— 

0 

Z 

5,  000 

5,000 

Wind,  knots 

-30 

30 

Azimuth,  deg 

— 

0 

Z 

4,  000 

4,000 

Wind,  knots 

-20 

20 

Azimuth,  deg 

— 

0 

The  final  surface  output  parameters  common  to  both  programs  were  in  agree- 
ment except  for  one  particular  aspect.  In  the  two-dimensional  program,  the  wind 
direction  reversal  takes  place  only  in  the  sense  that  the  wind  velocity  decreases 
from  positive  through  zero  to  negative  values.  Thus  the  program  retains  the  single 
vertical  plane  containing  the  cable,  a horizontal  plot  shows  a straight-line  projec- 
tion of  the  cable,  and  while  not  computed  or  printed,  a value  of  I equal  to  zero  is 
inferred.  Figure  6 illustrates  the  straight  line  obtained  in  an  X-Y  (H-I)  plot  of 
the  76.  006  output. 

Also  shown  is  the  77.  007  output  which  has  a final  value  of  1 equal  to  -7 19  ft., 

This  is  caused  by  the  rotation  of  the  wind  azimuth  from  180°  at  an  altitude  of 
10,  000  ft  to  0°  at  an  altitude  of  8000  ft.  Thus  the  cable  experiences  a side  force  in 
this  region  and  turns  in  azimuth.  As  can  be  noted  in  Figure  6,  most  of  the  bend  is 
completed  at  8000  ft  altitude  and  no  curvature  in  the  horizontal  plane  occurs  from 
that  point  to  the  surface.  This  effect  is  correct  for  the  conditions  specified  by  the 
particular  way  in  which  the  wind  was  entered.  If  a south  wind  diminishing  from 
60  knots  to  zero  followed  by  a north  wind  increasing  from  zero  to  15  knots  is  the 
true  condition,  then  additional  points  should  be  specified.  If  a single  point  at  8400  ft, 
0 knots  at  180°  were  specified  (this  is  the  zero  wind  velocity  intercept  between  10,  000 
and  8000)  no  bend  would  occur  above  8400  ft.  However,  a bend  would  then  take 
place  between  8400  ft  and  8000  ft  if  an  end  of  a cable  element  fell  within  that  area. 
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Therefore  to  reduce  the  probability  of  any  bend,  two  sets  of  entries  1 ft  apart  is 
suggested  for  these  cases  where  all  winds  are  in  the  same  or  opposite  direction. 
Foi  this  case,  they  would  be: 


z 

8400 

Wind 

0 

Azimuth 

180 

Z 

8399 

Wind 

0 

Azimuth 

0 

When  these  are  included  in  the  wind  table  and  a run  is  made,  there  is  exact  agree- 
ment with  the  two-dimensional  program;  that  is,  no  side  displacement  occurs. 
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Figure  6.  Plan  View  of  Cable,  Balloon  Origin, 
Test  1 


B.  Test  2—180*  Azimuth  Rotation 

In  Test  1 the  cable  moved  to  a minus  value  of  I when  a wind  azimuth  inter- 
polation between  180*  and  0"  was  called  for  between  10,000  ft  and  8000  ft.  This 
indicates  that  counterclockwise  rotation  of  wind  azimuth  will  occur  when  an  exact 
180®  turn  is  encountered.  To  show  this,  the  following  problem  was  established. 
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llulloon  Altitude, 

ft  MSL 

10,  000 

Surface  Altitude, 

ft  MSL 

4.  000 

Internal  Computations 

V.  (Yes) 

Cable  Diameter, 

in. 

0.  28 

Cable  Weight/ 1000  ft,  lb 

25.  0 

Element  Length  (K),  ft 

500.  0 

Balloon  Total  Force,  lb 

2000. 0 

Total  Force  Elev 

. Angle,  deg 

85.  0 

Wind  Field 

Test  2A 

Test  2B 

Test  2C 

Z 

10, 000 

10,  000 

10,  000 

Wind,  knots 

25 

25 

25 

Azimuth,  deg 

27  0 

270 

270 

Z 

9000 

9000 

9000 

Wind,  knots 

50 

50 

50 

Azimuth,  deg 

90 

89 

91 

Z 

4000 

4000 

4000 

Wind,  knots 

25 

25 

25 

Azlx'.-iuth,  deg 

90 

89 

91 

The  program  will  always  assume  a rotation  of  wind  azimuth  between  specified 
input  points  to  the  less  than  180*“  direction.  In  Figure  7,  the  H-I  plot  shows  that 
the  cable  responds  by  a move  to  the  right  (-1)  from  balloon  to  winch  indicating  a 
wind  from  the  right  quadrant  for  both  Teats  2A  and  Test  2C.  Therefore,  when  an 
exact  180°  reversal  is  presented  (2A),  the  data  indicates  that  the  wind  rotation 
will  be  counterclockwise  the  same  as  in  2C  where,  between  10,  000  and  9000  ft,  a 
179“  rotation  from  270“  through  180“  to  91“  is  known  to  occur.  In  Test  2B  the 
rotation  of  179“  would  proceed  from  270“  through  0“  to  89“  or  retain  winds  from 
the  left  quadrant  in  Figure  7.  In  this  way,  very  different  winch  locations  are 
indicated. 

In  summation,  this  demonstration  shows  the  desirability  of  never  specifying 
two  adjacent  wind  azimuths  exactly  180“  apart.  If  there  is  any  knowledge  that 
would  aid  in  better  defining  the  direction  of  rotation  an  intermediate  point  should 
be  part  of  the  wind  field  input. 
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Figure  7.  Plan  View  of  Cable 
Balloon  Origin,  X„-Y„  Axes, 
Test  2 ® ^ 


4.  1.  2 EFFECT  OF  SELECTED  VALUE  OF  ELEMENT 
LENGTH-TEST  3 

Unless  a value  for  the  length  of  the  cable  element,  K,  used  in  the  computations 
as  specified,  a value  equal  to  1 percent  of  the  height  or  difference  between  the 
balloon  and  surface  altitude  is  program- selected.  In  order  to  determine  a reason- 
able value  of  K such  that  the  number  of  program  interactions  and  running  time  be 
minimized,  a series  of  runs  was  made  with  K varying  between  0.  1 of  1 percent  to 
10  percent  of  the  height.  While  this  was  done  for  only  one  specific  problem— height, 
balloon  force,  cable  specifications,  and  wind  field— it  does  indicate  something 
about  the  desirable  order  of  magnitude  of  the  K value. 
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1’est  3 consisted  of  the  following  entries: 


Balloon  Altitude,  ft  MSL 
Surface  Altitude,  ft  MSL 
Internal  Computations 

Cable  Diameter,  in. 

Cable  Weight/ 1000  ft,  lb 
Element  Length  (K),  ft 

Balloon  Total  Force,  lb 
Total  Force  Elev.  Angle,  deg 


Wind  Field 

"Z  14.000 

Wind,  knots  25 

Azimuth,  deg  180 

Z 13,000 

Wind,  knots  25 

Azimuth,  deg  225 

Z 11,000 

Wind,  knots  35 

Azimuth,  deg  270 


14,000 
4,000 
V (Yea) 

0.  28 


25.0 

Varies,  see  below 

3,  200.0 

85.0 

Z 

10, 000 

Wind,  knots 

40 

Azimuth,  deg 

300 

Z 

8500 

Wind,  knots 

30 

Azimuth,  deg 

270 

Z 

4000 

Wind,  knots 

15 

Azimuth,  deg 

210 

K Values 


ft 

% of  Zg-Zg 


1000 

10 


500  400  200  100 
5 4 2 1 


50 
. 1 


10 
. 1 


Each  of  the  K values  was  used  in  a run  with  three  different  printouts.  Mag- 
netic program  cards  were  made  up  in  three  different  forms  as  follows: 

(a)  Print  of  final  output  at  the  surface  only, 

(b)  Print  of  Z,  H,  and  I at  end  of  each  cable  element  plus 

final  output  at  the  surface, 

(o)  Complete  print  of  program  output  as  written— all  parameters 
at  end  of  each  cable  element  plus  final  output  at  the  surface. 

The  time  required  to  run;  (a)  each  of  the  21  combinations  without  plotting 
(77.  007),  and  (b)  some  combinations  with  plotting  (77.  007P)  was  also  measured. 
Table  3 presents  the  values  of  several  computed  parameters.  The  values  for  a K 
value  of  10  ft  (0.  1 of  1 percent  of  Zg-Zg)  were  taken  as  having  zero  error  and  were 
used  as  a base  for  error  determination  for  other  K values.  These  errors,  shown 
in  Figure  8,  indicate  that  a K value  of  approximately  3 percent  of  Zg-Zg  would 
assure  less  than  1 percent  error  in  most  of  the  parameters.  The  lack  of  smooth- 
ness in  the  error  curves  can  be  attributed  to  the  manner  in  which  the  final  surface 
altitude  is  reached.  A function  of  the  K length  and  cable  elevation  angle,  the  sur- 
face intercept  can  show  somewhat  random  error.  This  in  turn  can  introduce  addi- 
tional smaller  errors  in  other  parameters. 

However,  the  general  trends  are  indicative  of  the  error  in  using  too  large  a 
K value  such  as  10  percent  of  Zg-Zg. 
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K , ft, 

10 

f,0 

100 

800 

400 

800 

1000 

thivf  ,Mt,. 

ft 

4000,081 

4000,884 

4000,418 

4014,808 

4014,181 

4046,488 

4041,461 

'I'anstoii 

It) 

8480,16% 

3440.404 

8880.641 

8081,681 

8488,484 

8416,011 

8414,866 

Ott^Uv  tAnqlti 

48,466 

48,464 

48,448 

48,466 

44,008 

44,040 

44,146 

C.Wsiuhl, 

in 

aA9,%8% 

848.400 

'148.800 

848,000 

888,000 

881 ,880 

880,000 

e,t,eiAqth 

ft 

10101,000 

10100,000 

10100,000 

10080,000 

lOOIO.OOO 

10040,000 

10000,000 

V.nt'atj 

Ih 

14,664 

14,414 

14,148 

11,448 

14,818 

16,418 

18,480 

n 

ft 

1014.414 

1014,080 

10,14.846 

1004,444 

481,444 

488.801 

411,064 

j 

ft 

884,011 

884,448 

884,448 

881,048 

114,886 

114,816 

101,618 

ft 

1111.441 

1104.444 

1104.184 

1848.811 

1888,088 

1848,444 

1810.486 

1 M 

18.484 

14.018 

14,110 

14,881 

14.818 

14,688 

40,644 

a, Kiev  snq 

1 

88,484 

88.414 

88,481 

88,884 

88,644 

88,480 

81 , Range 

ft 

10084,440 

10084,446 

10014.461 

10064,611 

10064,614 

10014, 4&4 

4484.648 

CaAniroufeh  49 

68.888 

81,081 

81,048 

81,188 

61.481 

81,610 

84,118 

Figure  W also  Indioatea  that  the  optioital  built  In  I pei'oent  Is  aeleotor  probably 
aaeurea  that  all  parametera  are  within  I pernent  of  true  valuea,  althoughi  aa  neat 
Indloated,  running  tlmoa  may  bo  eaoeaalve  for  oertaln  oomblnatlona  of  operation, 
Figure  9,  presents  the  program  running  times  aa  a function  of  the  eable  element 
length,  K.  If  a K value  of  S percent  of  were  aelet.eil  from  error  conaldera- 

tlona,  the  following  running  times  would  be  Indicated, 


(al 

No  Plot, 

Only  h'lnal  Print 

a min 

(b) 

No  Plot, 

Z,  H,  1 Intermediate  Alt 
Points  and  Final  Print 

a -3/4  min 

(o) 

No  Plot, 

Complete  Intermediate 

Points  and  Final  Print 

7 min 

(d) 

Plot, 

Only  Final  Print 

6 min 

(e) 

Plot, 

7-,,  H,  1 Intermediate  \lt 

Points  and  Final  Print 

a min 

<f) 

Plot, 

Complete  Intermediate 

Points  and  Final  Print 

9-3/ 4 min 

It  should  be  emphaslaod  that  these  errors  and  running  times  apply  to  the  par- 
ticular problem  used  In  this  model.  \V;  P y differing  problems,  such  as  more 
variations  in  the  wind  aelmutb  or  greater  wind  magnitude  producing  larger— or 
smaller  cable  elevation  angles  at  the  surface— will  have  a different  set  of  errors 
and  of  course,  running  times. 


4.  1.3  SOLUTION  CHECK  FOR  ALL  BALLOON  POINTING 
AZIMUTHS -LARGE  CABLE-TESTS  4 THROUGH  12 


The  cables  specified  In  previous  tests  have  small  diameters.  They  are  typical 
of  the  type  of  cables  selected  for  tethered  balloon  operation  with  a high  strength- 
to-welght  ratio  so  that  small  dlEimeters  are  possible  thereby  minimizing  drag  loads 
as  well  as  cable  weight.  With  small  drag  loading,  the  cable  is  less  sensitive  to 
winds  and  in  particular,  winds  from  the  side  will  not  produce  large  deflections 
easily  shown  on  an  X-Y  plot.  In  order  to  clearly  show  large  movement  of  the  cable 
and  illustrate  that  this  program  can  properly  handle  large  cable  turns,  a peculiar 
type  of  cable  will  be  specified  in  this  group  of  tests.  In  addition,  the  wind  field  will 
include  some  large  wind  speeds. 


A large  cable  diameter  of  1.  25  in.  together  with  winds  of  large  magnitudes  will 
produce  significant  drag.  A cable  weight  of  only  20  lb  per  1000  ft  insures  that  cable 


weight  will  not  predominate  in  the  calculations  of  cable  deflections.  (This  is  not  a 


real  cable. ) Other  input  data  include; 

Zg  10,000  ft  MSL 

Zg  4.  000  ft  MSL 

Internal  ti  (Yes) 

K 180  ft 

F.J,  3,200  1b 

0 85  deg 

and  the  ilowing  wind  profiles; 


' de 

Wind 

Speed 

Test 

Test 

A 

Test 

z i m 
Test 

u t h. 
Test 

deg 

Test 

Test 

Test 

Test 

i -VISL 

knots 

4 

5 

6 

7 

8 

9 

10 

11 

12 

10,  000 

25 

180 

180 

180 

270 

270 

0 

0 

90 

90 

8,  500 

40 

180 

225 

135 

315 

225 

45 

315 

135 

45 

7,  000 

50 

180 

270 

90 

0 

180 

90 

270 

180 

0 

6,  000 

60 

180 

300 

60 

30 

150 

120 

240 

210 

330 

4,  000 

20 

180 

315 

45 

45 

135 

135 

225 

225 

315 

Test  4 is  a two-dimensional  problem  since  the  complete  wind  profile  shows  wind 
from  the  soirth  (180°).  The  balloon  therefore  points  towards  180“  and  the  winch  lies 
south  of  the  balloon.  Hence  the  azimuth  from  the  winch  to  the  balloon  and  the  azi- 
muth of  the  cable  leaving  the  winch  should  be  0°  as  was  computed  (see  Table  4). 

The  same  problem  could  be  run  in  Program  7 6.006  since  it  is  a two-dimensional 
condition. 
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Table  4.  Output  Parameters —Tests  4 Through  12 


Test 

Surface 

Height 

Tension r 

Cable  El. 

Cable 

Cable 

Tot. Vert, 

Bain. El. 

Slant 

Distance 

No. 

Altitude 

Winch 

Angle 

Length 

Weight 

Drag 

Angle 

Range 

H 

ft 

ft 

lb 

deg 

ft 

lb 

lb 

dog 

ft 

ft 

4 

4C05.25 

5994.75 

3124.00 

28.89 

8046 

160.92 

1517,72 

52.20 

7586.43 

4649.39 

5-12 

4003.69 

5996.11 

3271.11 

35.93 

7290 

145.80 

1122.26 

60.85 

6865.74 

1476.21 

” 

Distance 

Distance 

Azimuth  to 

Cable  Out 

X 

y 

I 

L 

Balloon 

Azimuth 

+ East 

+ North 

ft 

ft 

deg 

deg 

ft 

ft 

4 

0 

4649.39 

0 

0 

0 

4649.39 

5 

3000.98 

3344.42 

63.81 

85.55 

3000.96 

1476.21 

6 

-3000,90 

3344.42 

296.19 

274  .45 

-3000.98 

1476.21 

7 

3000.98 

3344.42 

153.81 

175.55 

1476.21 

-3000.98 

8 

-3000.98 

33*4.42 

26.19 

4.45 

1476.21 

3000.98 

9 

3000,98 

3344.42 

243.81 

265.55 

-3000.98 

-1476.21 

10 

-3000.98 

3344.42 

116.19 

94.45 

3000.98 

-1476.21 

11 

3000.98 

3344 .42 

333.81 

355.55 

-1476.21 

3000.98 

12 

-3000.98 

3344.42 

206.19 

184 .45 

-1476.21 

-3000.98 

While  the  wind  speed  profiles  for  Tests  4 through  12  are  all  the  same,  the 
azimuth  profiles  of  Test  5 through  Test  12  indicate  winds  turning  in  azimuth  as 
the  altitude  decreases.  In  Tests  5 through  12  there  is  essentially  the  same  degree 
of  turning  in  order  to  represent  the  balloon  pointing  in  4 different  directions  with 
the  cable  subjected  to  winds  approaching  90“  from  either  side  of  the  cable.  The 
wind  azimuths  in  Tests  6,  8,  10,  and  12  are  mirror  images  of  those  in  Tests  5, 

7,  9,  and  11  respectively.  They  were  selected  in  order  to  show  positive  and 
negative  values  of  the  displacement,  I,  and  to  produce  a wide  range  of  cable  and 
balloon  azimuth  angles  and  X and  Y values  on  the  surface.  (Section  3.5,7  contains 
the  tape  output  for  Test  5 only. ) Table  4 contains  the  final  outpirt  values  to  Illustrate; 
(1)  parameters  such  as  tension,  cable  elevation  angle,  cable  length,  slant  range, 
etc, , are  identical  for  all  tests  as  they  should  be;  and  (2)  that  balloon  and  cable 
azimuth  angles  and  the  geographic -axis  displacements  X and  Y differ  as  expected. 
Figure  10  illustrates  the  eight  H-I  or  X^-Yg  plots  of  these  tests  to  show  that 
the  parameters —azimuth  to  the  balloon  from  the  winch  and  the  azimuth  of  the  cable 
leaving  the  winch— are  consistent,  and  therefore  properly  handled  through  the  0“ 
to  360®  boundary. 


Figure  10.  Plan  Views  of  Cables,  Balloon 
Origin,  Xg-Yg  Axes,  Tests  5 Through  12 

4.  1,  4 TRENDS  IN  CABLE  ROTATION  DUE  TO  WIND  AZIMUTH 
ROTATION- LARGE  CABLE-TESTS  13  THROUGH  18 

The  balloon  and  cable  specifications  and  the  magnitude  of  the  wind  used  in 
Tests  4 through  12  were  retained  in  Tests  13  through  18.  As  shown  in  Figure  11, 
the  azimuth  angles  of  the  wind  on  the  cable  were  varied;  (a)  in  Tests  13  through  16 
to  greater  amounts  of  clockwise  rotation  than  Test  5,  and  (b)  in  Tests  17  and  18 
through  less  severe  rotation  clockwise  then  counterclockwise. 

The  H-I  (Xt3-Y„)  plot  In  Figure  12  indicates  the  tightness  of  the  cable  turn 
r5  xD 

produced  by  very  severe  wind  azimuth  rotation.  When  Figure  12  is  transferred  in 
Figure  13  to  a common  winch  set  of  axes,  X^-Y^,  a clearer  picture  of  the  balloon 
movement  with  variations  in  wind  rotation  on  the  cable  is  possible.  Because  the 
balloon  in  these  particular  tests  is  pointing  exactly  south  (180“  azimuth),  the 
Xw'^w  axes  in  Figure  13  are  also  the  X-Y  or  geographical  axes  with  North  pointing 
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up  Iho  Y-axos.  It  can  be  seen  that  the  *[roatest  amount  of  wind  rotation  (Teat  10) 
places  the  balloon  closest  to  the  winoh  in  this  horiviontal  projection  of  the  cable 
a’oometry.  The  winoh/balloons  in  Testa  17  and  18  and  in  Tost  4,  the  two-dimensional 


case,  are  outside  the  boundaries  of  Figures  12  and  13,  The  length  of  the  cable  re- 


quirod  in  each  of  the  tests. 

as  shown  below,  decreases  with 

proximity  to  the  winch 

as  projected  in  the  horizontal  plane. 

Horizontal  Distance 

Test  Number 

Balloon  to  Winch 

Length  Cable 

4 

4649 

8046 

18 

4445 

7920 

17 

4285 

7812 

5 

3344 

7290 

13 

2425 

6786 

14 

1192 

6300 

15 

828 

6210 

16 

515 

6120 

Min  Possible  C. 

Length 

Zb-^s 

6000 

The  two-dimensional  case,  winds  all  from  the  same  direction  (Test  4),  produces 
the  greatest  displacement  and  cable  length. 


Figure  11.  Wind  Profiles,  Tests  4,  5,  and  13  Through  18 
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Figure  14  presents  the  variation  of  tension  and  elevation  angle  with  altitude  for 
some  of  these  tests.  Test  5 maintains  the  highest  level  of  tension  of  all  the  tests. 
In  this  case  the  tension  at  the  winch  is  greater  than  at  the  balloon.  Its  elevation 
angle  at  the  winch  is  only  36".  While  Test  16  produces  a very  vertically  aligned 
cable  (at  6000  ft  it  is  nearly  vertical)  the  tension  decreases  down  the  cable  to  a 
minimum  at  the  winch.  The  more  tiian  8000  ft  of  cable  in  Teat  4 enters  the  winch 
at  a very  low  29"  elevation  angle.  While  these  effects  are  all  exaggerated  cases 
due  to  the  very  lightweight  large-diameter  cable,  they  illustrate  the  necessity  to 
make  computations  over  a wide  range  of  possible  meteorological  conditions,  if 
precise  information  on  cable  behavior  is  wanted  for  many  different  types  of  flying 
days. 


Figure  14,  Cable  Tension  and  Elevation  Angle  Variation  with  Altitude, 
Tests  4,  5,  13,  and  16 


4.  1.  5 EXTREME  CABLE  ROTATION  WITH  BALLOON  IN 

ALL  QUADRANTS-LOW  ALTITUDE  CYCLES-TEST  19 

The  following  problem  is  Introduced  to  illustrate;  (a)  a corkscrew-like  cable 
geometry  which  extends  through  all  4 H-1  quadrants,  and  (b)  the  use  of  a lower  alti- 
tude cycles  for  static  evaluation  of  the  balloon  during  ascent  or  decent.  Again,  a 
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large  diameter  and  extremely  lightweight  cable  is  utili^ed  to  obtain  the  exaggeration 
desired  and  is  not  an  available  product. 


Max.  Balloon  Alt,  Z 
Surface  Alt,  Zg 

Cable  Diameter 
Weight  per  1000  ft 
Element  Length  K 
Balloon  Total  Force 


B 


10,  000  ft  MSL 
0 ft  MSL 

1.  0 

2.  0 in. 

10  lb 
250  ft 

3200  lb  (also  used  at  all  lower  alts) 


Angle 

of  Total  Force 

85° 

Z 

10, 000 

7000 

4000 

1000 

Wind 

25 

40 

40 

40 

AZ 

180 

45 

180 

315 

Z 

9000 

6000 

3000 

0 

Wind 

35 

40 

40 

40 

AZ 

27  0 

90 

225 

0 

Z 

8000 

5000 

2000 

Wind 

40 

40 

40 

AZ 

0 

135 

270 

Note  that  the  wind  magnitude  from  8000  ft  to  the  surface  is  40  knots  and  that  a con- 
stant Cq  = 1.  0 is  used  in  this  problem. 

Figure  15  istheH-I  (Xg-Yg)  plot  withthe  balloon  at  10,  000  ft  showing  that  the  cable 
extends  in  a spiral  through  all  4 H-I  quadrants.  Figure  16  shows  the  Z-H  and  Z-I 
vertical  views  of  the  cable  (at  a different  scale  than  Figure  15)  to  show  the  tightness 
of  the  spiral. 

Figure  15  was  then  converted  to  an  X^-Y^  plot,  Figure  17a,  Again,  because 
the  azimuth  of  the  wind  at  the  balloon  at  10,  000  ft  is  180°,  these  axes  can  be  con- 
sidered as  X-Y  with  North  up  the  paper.  The  cable  for  Zg  = 10,  000  ft  is  shown 
appearing  smaller  than  in  Figure  15  due  to  the  reduced  scale. 

The  X and  Y values  from  the  lower  balloon  altitude  cycle  runs  (Option  2)  were 
then  plotted  as  points  on  Figure  17a.  The  wind  azimuth  at  the  balloon  at  each  of  the 
1000  ft  levels  was  used  to  "aim"  the  small  balloons  (drawn  thereon)  into  the  wind. 
Each  is  annotated  with  balloon  altitude.  As  would  be  expected  with  the  wind  field 
established  for  this  test,  the  balloon  travels  through  a spiral  path  relative  to  the 
winch  during  ascent  from  the  surface  to  10,  000  ft. 

In  the  Interests  of  clarity  in  Figure  17,  the  cable  plan  views  are  drawn  only 
for  balloon  altitudes  of  10,000,  8000,  6000,  and  4000  ft.  In  spite  of  the  widely 
dispersed  balloon  locations,  the  cable  leaving  the  winch  remains  within  a 45°  to 
155°  range  of  azimuths  at  all  balloon  altitude.  Figure  17b. 

This  type  of  exaggerated  case  also  causes  the  cable  tension  to  increase  moving 
from  balloon  to  surface.  Figure  18a.  At  all  lower  balloon  altitudes,  the  winch 
tension  is  larger  than  the  starting  tension,  3200  lb.  The  largest  tension  is  at  the 
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q,  BALLOON  AT  10,000  ft,  b,  BALLOON  AT  INDICATED  ALTITUDE, 

CABLE  TENSIOtJ  WINCH  TENSION 


Figure  18.  Tenalon— Test  IS 


4.  1,  6 fHANGES  WHEN  LESS  EXTHEME  CABLES  OH  AZIMUTH 
HOTATIONS  AHE  INTRODUCED-TESTS  19  THROUGH  24 

The  problont  in  Test  19  used  u severely  rotating  wind  field  and  a large  liglit- 
woight  cable  to  illustrate  a spiral  rotation  of  the  cable.  Tests  20  through  24  retain 
thu  same  wind  magnitudes  us  in  Test  19  but  vary  other  parameters. 

Tost  20  uses  the  same  cable  but  is  a two-dimensional  case  with  the  winds  at  all 
levels  from  the  same  direction.  In  Figure  10,  the  effect  may  be  noted  as  large 
increases  in  downwind  displacement  and  in  cable  length  and  a very  low  cable  eleva- 
tion angle  at  the  wlncln 

Tests  23  and  21  are  three-dimensional  and  two-dimensional  repeats  of  the  same 
problem  as  Test  19  and  20  except  with  a smaller  diameter  cable  (0.  5 in.  ) to  reduce 
the  drag  components.  While  winch  tensions  show  a small  decrease  ewer  Tests  19 
and  20,  tile  elevation  angles  are  increased  and  tile  cable  length,  (in  the  two-dimen- 
sional case)  is  greatly  reduced. 
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Figure  19.  Vertical  Views  of  Cable,  Tests  19  Through  24 


Tests  24  and  22  are  three-dimensional  and  two-dimensional  repeats  of  the 
same  problems  except  that  the  cable  is  brought  into  a completely  realistic  specifica- 
tion by  an  increase  in  weight  to  100  lb  per  1000  ft  with  a diameter  of  0.  5 inches.  In 
the  comparison  of  two-dimensional  cases,  the  cable  length  and  down  range  displace- 
ment increase,  but  the  winch  tension  is  significantly  reduced.  In  tlie  comparison  of 
three-dimensional  cases,  the  winch  tension  is  similarly  reduced  without  much 
change  in  cable  elevation  angle. 

Tests  23  imd  24  both  exhibit  a different  form  of  cable  rotation  in  the  plan  view 
than  Test  19.  The  effect  of  the  smaller  drag  producing  cable  is  similar  in  both 
Tests  23  and  24,  therefore  Figure  20  Includes  results  for  only  Test  24.  In  the  H-1 
plot— left  side— the  cable  becomes  vertical  at  some  altitude  and  produces  the  sharp 
discontinuity  shown  thereon.  This  is  the  special  case,  discussed  in  Section  3.  2.  8b 
which  required  special  computational  handling  within  the  programs.  On  the  right 
side  of  h'igure  20,  a portion  of  the  vertical-plane  plots  of  Z-H  and  Z-I  show  that  the 
cable  goes  vertical  at  an  altitude  of  approximately  6600  ft.  These  comparisons 
indicate  that  a practical  tether  cable  will  rarely  produce  the  spiral  cable  configura- 
tion in  a wind  field  showing  continuous  rotation  with  changes  in  altitude.  More  likely 
is  the  form  shown  in  Figure  20  with,  in  addition,  a decrease  in  cable  tension  as 
one  moves  down  the  cable  to  the  surface,  as  with  Test  24. 
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Kigui'o  20,  Plan  and  Vortical  Views  of  Cnblo,  Tost  24 


4,2  Prublem*  fur  Totting  Program  77.007U  Opoiatlon 

When  Tost  24  was  completed  in  a Program  77,  007  run,  Option  3 was  selected 
to  call, in  Program  77,  007 B.  This  ostublished  a fixed  cable  length  of  10,  073  ft 
which  was  found  in  the  Test  24  solution  (Col,  1 below).  The  wind  profile  for  Teat 
22,  a two-dimensional  case— wind  ail  from  same  direction  — with  wind  magi\ltudes 
the  same  as  Tests  10  through  24,  was  introduced  ns  an  input  into  Program  77. 007  B. 
This  test  will  therefore  be  designated  as  Test  24.  22, 

Some  of  the  output  is  shown  in  Col.  2 below.  It  indicates  that  the  balloon  has 
descended  from  an  altitude  of  10,000  ft  (Test  24)  to  8980  ft  (Test  24.  22)  in  the 
presence  of  the  two-dimensional  wind.  To  indicate  consistency  within  the  programs, 
the  altitude  of  0980  ft  was  next  entered  in  Program  77.  007  and  a solution  made  again 
with  same  cable  and  winds  used  with  77.  007  B.  The  surface  output.  Col.  3,  exactly 
duplicates  the  Program  77,007B  output,  If  a printout  and/or  plot  of  cable  param- 
eters vs  altitude  were  required,  this  is  the  procedure  to  follow  after  a 77,  007  B 
solution. 
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(1) 

(2) 

(3) 

Test  24 

77. 007 

Test  24,  22 

77. 007 B 

Test  22.  22 
77.  007 

Program  No.  : 


Input  Data 


Balloon  Alt. 

10, 000 

— 

8,  980 

Winds 

3-dim. 

2-dim. 

2-dim. 

Cable  Length 

— 

10, 075 

— 

F 'j' 

3.  200 

3,  200 

3,  200 

e 

85 

85 

85 

Output  Data 

Balloon  Alt. 

— 

8,  980 

— 

Cable  Length 

10,  075 

10.  075 

10, 075 

Winch  Tension 

2,  219 

2,  322 

2.  322 

Cable  Ele-v.  Angle 

82.  0 

42.  8 

42,  8 

Azim.  to  Bain. 

352,  0 

0 

0 

Cable  Azimuth 

48.  8 

0 

0 

X 

-139 

0 

0 

Y 

990 

4,  054 

4,  054 

4.3  Practical  Problems— High  Altitude  Tethered  Balloon 

Since  the  developments  of  these  programs  were  predicated  on  their  need  in 
solving  advanced  tethered  balloon  problems,  some  illustrative  cases  are  pertinent 
here.  As  described  in  Reference  1,  there  is  an  AFGL  plan  to  tether  a balloon  at 
an  altitude  of  20  km  (65,616  ft  MSL),  Two-dimensional  cable  problems  and  their 
solutions  were  presented  in  Reference  1,  Section  5,  utilizing  Program  No.  76.  006. 

4.  3.  1 FIXED  BALLOON  ALTITUDE 

To  illustrate  the  effect  of  winds  having  directions  that  vary  with  altitude  on  the 
cable  design  and  system  performance,  one  of  the  aforementioned  two-dimensional 
problems  will  be  repeated  and  expanded  here  into  three  dimensions.  The  basic 
problem  using  the  design  wind  magnitude  profile,  the  internal  drag  coefficient 
computations,  and  the  same  cable  as  shown  in  Reference  1,  Table  3,  Line  3 was 
selected  as  a typical  high-altitude  problem.  Test  25,  with  its  input  shown  below, 
is  a repeat  of  that  problem  using  wind  azimuths  at  all  altitudes  equal  to  27  0°  to 
provide  the  comparable  two-dimensional  case  in  Program  77.  007. 
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Input  - Teat  25 


a 

65,616  ft 

Z 

= 

65,616  ft 

Z 

= 

36, 000 

s 

4400  ft 

Wd 

s 

25  knots 

Wd 

= 

26 

Olam 

= 

0.  3 In. 

AZ 

= 

270  deg 

AZ 

= 

270 

Wt 

25.  0 lb/ 100  ft 

Z 

_ 

52,500 

Z 

29, 500 

K 

= 

1000  ft 

Wd 

15 

Wd 

— 

18 

Ft 

= 

3073  lb 

AZ 

270 

AZ 

= 

270 

e 

= 

81.6  deg 

Z 

42,  500 

Z 

S 

5500 

Wd 

= 

27.  5 

Wd 

= 

9 

AZ 

= 

27  0 

AZ 

= 

270 

Z 

= 

4400 

Wd 

= 

8 

AZ 

= 

270 

The  output  parameters,  shown  in  the  first  line  of  Table  5A  agree  with  those  of 
Reference  1.  Tests  26  through  30  were  then  run  with  the  same  input  as  Test  25 
except  for  wind  azimuth.  The  azimuths.  Figure  21,  were  widely  varied  from  test 
to  test:  those  used  in  Test  29  represent  a greater  than  360°  continuous  azimuth 
rotation  between  the  balloon  and  the  surface. 

WIND  SPEED  (knot*) 


(•'igurc  21.  Wind  Profiles,  Tests  25  Through  30 

In  Tabic  5A  it  may  be  notbd  that  the  widely  varying  azimuth  conditions  intro- 
duced in  Tests  26  through  30  produced  little  change  in  the  tension  at  the  winch  but 
increased  the  cable  elevation  angle  at  the  winch  seemingly  in  proportion  to  the 
"severity”  of  the  azimuth  rotation.  Cable  length  was  decreased  by  small  amounts 
but  the  balloon  decreased  Its  horizontal  displacement  from  the  winch  by  as  much 
IS  50  percent.  I'or  a given  fixed  profile  of  wind  magnitutle,  the  changes  in  wind 
fllrection  act  principally  to  Improve  (Increase)  cable  elevation  angle  at  the  winch 
aofl  to  improve  (decrease)  balloon  displacement. 
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Table  5.  High-Altitude  Tether-Cable  Test  Data— Tests  25  Through  30 


10'! 


32,3  298  1767  51.5  " 135.8  65.4  62110  145.4  18057  -18557 


l''or  the  condition  of  fixed  balloon  altitude— cable  len^h  varied  to  maintain  the 
height— and  fixed  balloon  wind  magnitude*  two  conclusions  seem  to  be  in  order, 

(a)  If  the  cable  "reaches"  the  surface  when  any  two-dimensional  wind  profile 
I'eonstant  aaimuth)  typo  of  problem  is  solved,  the  tension  at  the  winch  will  never 
be)  more  tiiun  a few  percent  higher  than  that  computed  for  a zero-wind  case.  In 
Koforenoa  1,  for  the  same  balloon  and  cable,  various  wind  magnitude  profiles  on 
the  cable  produced  tlie  same  tension  at  the  winch  in  spite  of  widely  varying  cable 
lengths,  total  cable  weights,  and  drag  forces.  Therefore,  the  tension  calculated 
with  the  simple  no- wind  (on  the  cable)  relationship 

'^W  " Cable “ F.p  -(Length  X Cable  Density) 


or 


“ F.J,  - (Zg-Zg)  X Cable  Density 


or 


n 3073  - (65,  618  - 4400)  X (25/1000)  = 1542.  6 lb 

is  within  1 lb  of  the  tension  computed  in  Tost  25. 

(b)  When  in  any  two-dimensional  wind  profile  type  of  problem  the  total  vertical 
drag  is  found  to  be  less  than  the  total  cable  weight,  a change  to  any  azimuth  profile 
(three-dimensional  with  same  wind  magnitude)  will  cause  little  significant  change  in 
the  winch  tension  provided  the  cable  "reaches"  the  surface.  There  are  even  some 
combinations  of  unusual  cables  and  heights  investigated  (see  Tests  4 and  5,  Table  4) 
where  the  vertical  drag  may  exceed  the  cable  weight  without  the  wind  azimuth 
variations  causing  excessive  adverse  tension  changes. 

Increases  in  the  wind  magnitude  at  the  balloon  can  of  course,  change  the  balloon 
total  force  and  therefore,  the  winch  tension  and  modify  the  above  statements;  as 
will  be  shown  in  Section  4.  3.  2b. 

4.  3.  2 FIXED  CABLE  LENGTH 
(a)  Effect  of  Wind  Azimuth  Change 

If  the  cable  length,  62,  600  ft  found  for  the  conditions  specified  in  Test  30  above 
is  held  fixed,  the  winds  (varying  only  In  azimuth)  from  some  of  the  previous  tests 
may  be  introduced  through  77,  007B  to  determine  resulting  new  balloon  altitudes 
and  cable  conditions.  Test  30.  29  in  Table  5B,  for  example,  is  a combination  of  the 
Test  30  cable  length  with  Test  29  winds.  The  balloon  total  force  and  angle  used  in 
these  tests  are  held  fixed  at  3073  and  81.6.  This  is  admittedly  imprecise  but 
sufficiently  close  enough  to  illustrate  trends  in  the  system  behavior  that  are  due 
principally  to  changes  in  the  wind  azimuths. 
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Table  5B  indicates  the  results  of  such  runs.  As  expected,  the  balloon  is  moved 
down  in  altitude  from  the  original  65,  616  ft  when  the  degree  of  rotation  of  the  wind 
between  the  balloon  and  the  surface  is  less  severe  than  the  "base"  wind  (Test  30). 
For  the  two-dimensional  v/lnd  case  (Test  30.  25),  the  balloon  drops  2300  ft  in  alti- 
tude. Its  downrange  displacement  nearly  doubles;  10,467  to  19,478  ft.  The  cable 
elevation  angle  at  the  winch  is  reduced  from  70  to  56°  with  a tension  increase  of 
59  lb.  When  a greater  amount  of  rotation  is  introduced,  the  balloon  rises  in  altitude 
as  shown  in  Test  30.  29;  a constant  25  knot  wind  magnitude  was  held  above  the  start- 
ing altitude  of  65,  616  ft. 

(b)  Effect  of  Wind  Magnitude  Change 

Here  the  cable  length,  62,  600  ft,  found  with  Test  30  conditions  is  again  held 
fixed.  Wind  magnitude  changes  were  made  as  follows: 


Z,  ft  Wind.  knots 


Original 

Teats  25 
Through  39 

Test 
30.  25 
30.  30 

66, 000 

25 

35.  36 

65, 616 

25 

35.  36 

52,500 

15 

26.72 

42,  500 

27.  5 

27,  5 

These  increases  at  52,500  ft  and  upward.  Figure  21,  represent  a doubling  of  the 

dynamic  pressure  and  resulting  drag  of  the  balloon  in  comparison  to  the  original 

values  at  65,  616  ft.  This  change  will  also  affect  parts  of  the  cable  at  these  altitudes. 

>1< 

The  wind  magnitudes  below  52,  500  ft  and  all  azimuth  values  were  left  unchanged. 

In  Table  5C,  Test  30.  25  uses  the  wind  magnitudes  shown  above  taken  with  the 
azimuths  of  Test  25  or  two-dimensional  case.  In  this  test,  the  balloon  moves  from 
63,  316  down  to  58,  260  ft.  In  Test  30.  30  the  balloon  level  changes  from  65,  616  to 
60,  920  ft.  In  both  cases  the  loss  in  altitude,  between  4700  to  5000  ft,  is  due  prin- 
cipally to  balloon  drag.  Both  of  these  tests  indicate  winch-tension  increases  of 
several  hundred  pounds  coupled  with  additional  decreases  in  cable  elevation  angle. 


In  this  exercise  the  balloon  total  force,  F-p,  and  its  angle,  9 , were  changed  to 
3170  lb  and  73.  5°  by  the  following  rationale.  With  the  assumption  of  a natural- 
shape  balloon  having  negligible  aerodynamic  lift,  when  the  dynamic  pressure  is 
doubled,  only  the  doubling  of  drag  was  considered.  Changes  of  size  and  shape  with 
altitude  were  considered  negligible. 

Orig.  Condition  New  Condition 


F.J,  = 3073  lb 


9 = 81.6° 
Drag,  D 


Frj,  = 3170  lb 
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9 = 73. 5° 

— ► D'  = 2D 


It  is  evident  that  the  increase  in  the  balloon  drag  (estimated  at  approximately  450 
lb)  is  felt  at  the  winch.  Changes  in  cable  drag  at  the  affected  upper  altitudes  also 
add  to  the  winch  tension  changes.  Increases  of  more  than  100  percent  in  the  down 
range  displacement  may  also  be  noted. 

The  results  above  indicate  that  the  cable  can  be  sized  principally  by  considera- 
tions of  basic  balloon  forces.  The  height,  and  cable  density  in  addition  will  in  most 
cases  determine  winch  tension  within  reasonable  working  tolerances.  Cable  eleva- 
tion angle  at  the  winch  becomes  a good  sensor  of  the  collective  effects  of  wind 
changes  on  the  cable  alone. 

The  test  cases  are  shown  here  to  illustrate  the  ease  of  using  the  programs  to 
solve  a typical  tethered-balloon  project's  design  and  flight  problems,  A great  many 
more  points  in  a matrix  of  variables  must  be  evaluated  to  cover  all  possible  condi- 
tions that  might  be  encountered  during,  for  example,  a flight  of  tworweek's  duration. 
Such  a parametric  study  for  the  high-altitude  tethered  balloon  would  also  have  to 
Include  the  effects  of  a changing  balloon  shape  and  size  during  ascent.  Unlike  the 
conventional  tethered-balloon  design  having  a ballonet,  this  balloon  is  a natural- 
shape  type  with  factory  installed  reefing  points  for  confining  the  excess  material  at 
altitudes  below  the  maximum.  As  the  balloon  rises,  one  reefing  point  at  a time  is 
released  to  permit  the  gas  to  expand  and  maintain  a non-flacid  shape.  Thus  a series 
of  values  of  balloon  total-force  F.^.,  and  its  angle,  6 , must  be  evaluated  for  many 
altitudes  to  provide  inputs  for  ascent  studies. 
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Appendix  A 

Symbols  and  Definitions 


aQ.  aj 

A 

AZ 
AZ 


B 


AZ, 


C 

C 

D 


D 


H 


D, 


g 

h 

H 

i 

I 

j 

J 

K 


Constants  used  in  Atmospheric  Density  Equation 

2 

Frontal  Area  of  Cable  Element,  ft 
Azimuth  of  the  Wind,  deg 

Azimuth  of  the  Balloon  from  the  Winch  Position,  deg 

Azimuth  of  the  Cable  Leaving  the  Winch,  deg 

Counter  in  77.  007 B 

Drag  Coefficient  of  a Cylinder 

Total  Wind  Drag  on  Cable  Element,  lb 

Horizontal  Component  of  D,j,  in  Vertical  Plane 
of  Element,  lb 

Horizontal  Component  of  D_  Perpendicular  to 
Vertical  Plane  of  Element, '^Ib 

Vertical  Component  of  D^,  lb 
Total  Balloon  Force,  lb 

Horizontal  Projection  of  Element  Length,  K,  ft 
Horizontal  Distance  along  Yg  axis,  ft 
Sum  of  Horizontal  Distances,  h,  ft 
Horizontal  Distance  along  Xg  axis,  ft 
Sum  of  Horizontal  Distances,  i,  ft 
Vertical  Projection  of  Element  Length,  K,  ft 
•Sum  of  Vertical  Distances,  j,  or  Zg  - Zg,  ft 
Incremental  Cable  Element  Length,  ft 
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L Horizontal  Distance  from  Winch  to  Balloon,  ft 

M Tension  Tick- Mark  Interval,  lb  (Plot) 

N Spacial  Tick-Mark  Interval,  ft  (Plot) 

nQD  Recall  Code  Number  in  Cj^  Computations 

NW  Recall  Code  Number  in  Wind  Computations 

ON  Option  Code  Number 

P Rounding  Factor  (Plot) 

2 2 

q Dynamic  Pressure,  l/2pV  , Ib/ft 

r Repeater  Code  Number  in  Optional  Lower  Altitude 

Runs  in  77.  007  and  77.  007P  or  in  77.  007B 

R Reynolds  Number 

T Tension,  lb 

V Wind  Velocity,  fps  ' 

Component  of  Wind  Velocity  In  line  with  Element,  fps 

Vj^  Component  of  Wind  Velocity  Normal  to  Element,  fps 

Vp  Horizontal  Component  of  Wind  Velocity  in  Vertical 

Plane  of  Element,  fps 

Vo  Horizontal  Component  of  Wind  Velocity  Perpendicular 

to  The  Vertical  Plane  of  the  Element,  fps 

W Weight  of  Cable  Element,  lb 

Wd  Wind  Velocity,  knots 

Xg  X-Axis  Centered  at  Balloon,  aligned  90°  with  Centerline 

of  Balloon  and  Positive  to  Right  of  Balloon  (View 
from  Above) 

Xy,  X-Axis  Centered  at  Winch,  parallel  with  Xg  axis  but 

Positive  Opposite  to  Xg 

X X-Axis  Centered  at  Winch,  aligned  East-West, 

Positive  East 

Yg  Y-Axis  Centered  at  Balloon,  aligned  with  Centerline 

of  Balloon  and  Positive  Forward  of  the  Balloon 

Y^  Y-Axis  Centered  at  Winch,  parallel  with  Y„  axis,  and 

Positive  Opposite  to  Yg 

Y Y-Axis  Centered  at  Winch,  aligned  North-South, 

Positive  North 

Z Altitude,  ft  MSL  and  Vertical  Axis 

Z„  Balloon  Altitude,  ft  MSL 

Id  ' 

Zg  Surface  Altitude,  ft  MSL 
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Greek  Symbols 


a Relative  Wind  Angle  to  V'ertical  Element  Plane,  deg 

^ Angle  of  Rotation  between  Adjacent  Vertical 

Element  Planes,  deg 

/3  „ Azimuth  of  the  Wind  at  the  Balloon  or  the  Azimuth 

to  which  the  Balloon  Points,  deg 

y Angle  in  the  X-Y  Plane  between  the  X-Axis  and  the 

Straight  Line  from  Winch  to  Balloon,  deg 

6 Angle  in  the  X-Y  Plane  between  the  X^-Axis  and  the 

Straight  Line  from  Winch  to  Balloon,  '^deg 

€ Elevation  Angle  of  Balloon  from  the  Winch,  deg 

6 Elevation  of  Element  or  Tension  Vector  above  the 

Horizontal,  deg 

(i  Atmospheric  Coefficient  of  Viscosity,  Ib/ft-sec 

3 3 

p Atmospheric  Density,  slugs /ft  (For  R,  Ib/ft  ) 

0 Angle  between  the  V and  Vectors,  deg 

■n  To  Call  Special  Operation— Computation  of  C^j 

A Increment 

T,  Sum 
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